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I. The Distribution of the Products of Heavy Ion Reactions with 
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By G. A. Cuackett, K. F. Cuackett and J. H. FREMLIN 
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[Received August 21, 1954] 


ABSTRACT 

The distribution in depth of the radioactive products of heavy ion 
bombardments of aluminium has been studied by stacked foil techniques. 
The bombardments were mainly done with 14N** ions of energies up to 
125 mev but a few were done with 160%+ ions. Some products were 
estimated from gross decay curves and some by radiochemical analysis. 
The results show that the products fall into two groups: *4Na together 
with heavier nuclei on the one hand, and 18F together with lighter nuclei 
on the other. In discussion it is shown that this can readily be explained 
if most products are derived from interaction of parts only of the nuclei 
concerned, in a way analogous to stripping reactions, rather than from 
breakdown of a compound nucleus. 


§1. INTRODUCTION 
In a previous paper (Chackett and Fremlin 1954) some experimental work 
was described on the production of N when aluminium was bombarded 
by 14N ions accelerated in the Birmingham cyclotron. It was shown that 
the greatest 4N activity occurred 60 microns from the surface of a thick 
target, indicating that the '*N nuclei produced retained a large part of the 
initial forward momentum of the 14N ions. This result was consistent 
with the formation of 1#N by loss of a neutron to the ?‘Al in a relatively 


weak interaction. 


* Communicated by the Authors, 
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It was believed that a similar investigation of the distribution in depth 
in the target of other reaction products would give useful information as to 
their mode of formation and in this paper we present the main results of 
such an investigation. 


§2. GENERAL MnrHop oF EXPERIMENT 


In all experiments, stacks of aluminium foil were bombarded by the 
internal beam of the cyclotron at the radius at which the maximum beam 
of ions with energies over about 50 Mev was obtained. For 14N ions this 
radius was 63:5 cm and for 160 it was 60cm. The target used, and the 
method of carrying out bombardments, have already been described 
(Chackett, Fremlin and Walker 1954, Chackett and Fremlin 1954). 
Since the earlier work the internal beams have been considerably improved, 
and measurable activities of the main products were observed down to a 
depth of 120 in the target, with a beam current of about 0-03 wa. 

Concurrently with the present work, the cyclotron was often engaged in 
the production of very large activities by deuteron bombardment. 
Consequently there was a serious risk of irrelevant activities being picked 
up by the stack of foils and, to prevent contamination, a 5y thick 
aluminium foil was always used to cover the main stack. This ‘ cover 
foil’ was discarded when there was any risk of confusion. 

The foils were separated immediately after bombardment, some being 
used for chemical separation of the isotope to be investigated, the activity 
of the rest being determined as a measure of the total effective bombard- 
ment. The strength of the bombardment could not be derived accurately 
from a current measurement as the total thickness of the stacks varied. 
Also small changes in the energy spectrum of the bombarding ions would 
give large relative changes in residual current traversing the whole stack. 
Relative bombardment intensities were therefore determined by measuring 
the **P activity of the otherwise unused foils after all other observable 
activities had decayed, and comparing these with the activity given by a 
standard bombardment. Comparisons of the *°P activity were nearly 
always accurate to better than 25°% and sometimes better than 10%. 
In all cases activities were observed with end-window Geiger counters, 
usually with a 7 mg/cm? window (G.M. 4). 


§3. Enercy DistrRIBuTION OF THE BEAMS 


Owing to the complex mechanism of acceleration of multiply-charged 
heavy ions (Walker, Fremlin, Link and Stephens 1954) the energy of the 
internal beam is not well defined as it is for deuterons. 

This makes interpretation of the results of the experiment to be described 
much more difficult but there is at present no way of avoiding the difficulty. 
An extracted beam has been obtained, but since it has not yet been possible 


to ‘shim’ the cyclotron appropriately for heavy ions, it is too small for 
radioactivation. 


Products of Heavy Ion Reactions with Aluminium 3 


The distribution of energy in the 14N beam at a radius of 63-5 cm under 
normal bombardment conditions was obtained as follows. The usual 
foil stack on the internal target was replaced by a thin platinum foil 
(2: mg/cm?) and 12 mm behind this, on the usual current collector, was 
placed an Ilford C2 nuclear emulsion. This was so arranged as to receive, 
at a small angle to its surface, '4N ions scattered elastically from the plat- 
inum foil together with any reaction products (mainly protons). The 
**N tracks in the developed plate were readily distinguishable by their 
appearance, and a range-histogram constructed for about 600 tracks at a 
scattering angle of about 40°. The range-energy relation is not accurately 
known but can be estimated well enough* for the present purpose and 
hence an energy-histogram can be obtained. For deflection at this angle 
from platinum, the Rutherford scattering law should be obeyed up to 
energies well beyond the potential barrier, viz. 80 Mev (lab. system) ; 
up to about 100 mev therefore the scattering probability for the given 
angle must vary as 1/H? where FH is the energy of the nitrogen nucleus. 
The numbers obtained in the energy-histogram must then be multiplied 
by ? to give the distribution of energy in the original beam. The 
appropriate small allowance for loss of energy in the platinum has also 
been made, giving finally the range-distribution shown in fig. 1. 


Fig. 1 
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Range distribution of 14N®* ions in aluminium. 


At this radius the highest energy to be expected is in the region of 
125 mev, but it can be seen that there is a maximum in the distribution at 
a little less than half this value. This is consistent with the measurements 
made here earlier with an internal crystal detector (Walker et al. 1954) 
although little of the latter work was done under the same operating 
conditions of the cyclotron. (In particular we have used somewhat 
higher gas pressures in the vacuum tank.) In considering the further 
results described below, it must of course be remembered that only 2/3 of 
the particle energy will be ‘ available ’ for reactions with aluminium. 


CONN i. 
* We are most grateful to Dr. D. Walker for the range-energy relation used. 


4 G. A. Chackett et al. on the Distribution of the 


$4. ResuLTS OF BOMBARDMENTS 


The distribution of yield of each of the reaction products from ™N 
bombardments which could be examined is shown in figs. 2 (a) and 2 (6d). 
The products in fig. 2 (a) (‘ heavy ’ products) comprise those of mass 24 to 
38, those in fig. 2 (b) (‘ light ’ products) of mass 11 to 18. 

In all cases the yields are given in terms of number of detectable atoms 
of each isotope actually coming to rest at in each 10 of aluminium at 
the depth concerned. The yields are normalized to a fixed yield of #*P 
found in the first 80. of the target below the 5 cover foil. Some of the 
results, namely those from 22P, ?8Al and the rather sketchy results for 1°O 
(the production of which in this way is reported here for the first time) were 
obtained entirely from analysis of the decay curves of different foils, 


Fig. 2 (a) 
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Distribution of ‘heavy’ products of 14N&+ bombardment of aluminium 
‘Random errors were significant only in the case of 24Na. The accuracy 
of normalization is usually better than 25%. 


using an elementary charge-discriminator when necessary. The results 
for *°S depend on differential counting using an end-window counter 
(window thickness 2 mg/em?) and a 47-geometry windowless counter. 
The majority, however, were obtained by chemical separation of the 
isotopes concerned. In all cases they are averaged from two or more 
experiments, the *°P and *4Cl curves being obtained by averaging results 
from a dozen or more experiments. 
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Most of the chemical methods used have already been described 
(Chackett, Fremlin and Walker 1954, Chackett and Fremlin 1954). 
In one or two cases, however, some improvements to the earlier methods 
_have been made. Thus, taking the elements in order of increasing Z, 
the following modifications were introduced : 

(i) Carbon (11C). The bombarded foils were dissolved in sodium 
hydroxide solution in a small reaction vessel through which passed a slow 
stream of air saturated at room temperature with ethanol vapour. The 
gas stream was roughly dried by bubbling through concentrated sulphuric 
acid and was then passed over copper oxide at 600°c. In this way the 
MC originally present in the foils was converted in the gas stream to 4CO, 
in the presence of a suitable quantity of CO, derived from the ethanol. 


Fig. 2 (b) 
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Distribution of ‘light ’ products of 14N¢+ bombardment of aluminium. N.b. 


(i) The curve for *2P is included for comparison. : 
: ; this 1 »e detecte 
il t oints for 15O are included ; this isotope could not be c 
ue ee Peicinty in the range 0-50 » and was not looked for beyond 70 p. 
The errors of normalization in this isotope are large. 


The gas was then bubbled through two centrifuge cones containing sodium 
hydroxide solution, which absorbed all the carbon dioxide. These cones 
were then detached from the apparatus and barium carbonate precipitated 
by addition of barium hydroxide. The barium carbonate was washed, 
mounted, and counted under an end window G.M. counters in the usu alway. 
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It carried the "C activity practically quantitatively and completely free 
from all other activities. 

(ii) Fluorine (18F). The method used previously for this isotope was to 
dissolve the aluminium foils in sodium hydroxide solution containing 
suitable carriers, then to remove the aluminium as the hydroxide by 
adjustment of the pH, and finally to precipitate the required activity, 
with carrier, as calcium fluoride. This method is not very satisfactory as 
much of the activity is lost on the aluminium hydroxide. It was found 
preferable to add finely divided silica to the solution (with carriers), acidify 
with perchloric acid, and distill off the fluorine as silicon tetrafluoride. 
The latter was trapped in cold water, the silicon removed by precipitation 
as silica by boiling with ammonium carbonate and the fluorine finally 
isolated as calcium fluoride or as lead chlorofluoride. 

(iii) Sodium (?4Na). The foils were dissolved in Toarceniten acid 
containing suitable carriers and sodium chloride was precipitated by 
addition of butanol saturated with hydrogen chloride. It was necessary 
to decontaminate the sodium chloride from other activities by repeated 
solution in a little hot water and precipitation by concentrated hydro- 
chloric acid at 0°c. At least six recrystallizations were performed. 

The activities induced by the bombardments were sufficient for a single 
operator to carry out chemical separations on six foils simultaneously for 
each element investigated having a half-life greater than half an hour, 
In the case of ?’Mg, 9°P and *8K it was rarely possible to process more than 
two samples simultaneously ; in the case of 4C the assistance of Mr. G. 
Clewes and Miss B. Hillman made it possible to make four separations 
from a single bombardment using the same apparatus for all of them. 


§ 5. Discussion 
5.1. Relative Total Yields of Different Isotopes 

We are now in a position to extend and amend the yield figures for the 
nitrogen reactions given in our earlier paper (Chackett et al. 1954), to 
take account of total yield rather than yield in the upper foils of the target 
only. ‘The figures are given in table 1 and are believed to be correct to 
25% for our standard conditions. Some particular ratios, such as 
30P/82P or 88K /94Cl are much better than this but there is no advantage at 
present in distinguishing these in view of the arbitrary nature of the 
conditions determining the original 14N energy-spectrum. 

The figures given are little changed apart from ®2P, 2°P and 85S which are 
all higher than had been thought before, because of the relatively greater 
yield of these isotopes in the lower foils and in the case of 35S, because of 
much more accurate counting. The preponderance of the neutron-excess 
isotopes is now even more striking. To explain this feature of the relative 
yields, we suggested that in many cases only part of the bombarding ion was 
captured by the target nucleus. The !4N nucleus could then be regarded 
as a group of quasi-independent particles any or all of which could be 
involved in reactions. On this ‘ buckshot ’ hypothesis, the distribution of 
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yields could be readily understood since a neutron would easily be captured, 
an alpha particle less easily and a single proton least easily owing to the 
considerable Coulomb barrier. 

A feature of importance in this process is that particles would be 
captured from a bound state to another bound state of very similar energy. 
There would therefore, be little excitation of the compound nucleus 
produced, apart from the energy acquired from the kinetic energy of the 
original **N. Hence it was shown that for the probable range of energies 
involved, there would usually be not more than one secondary heavy 
particle ejected for initial compound nuclei up to perhaps 3°Cl. 


Table 1. Relative Yields (in numbers of atoms) of Products of 14N 
Bombardments of Aluminium 


Isotope Present work Previous estimate 
(Chackett et al. 1954) 
el 500 300 
BeS 150 25 
28A] 140 140 
$401 100 (standard) 100 
Se KK 75 75 
Ane 60 . 25 
BON 20 g 
SiN 14 15 
24Na Il 12 
a2 0) 5 a 
23h 20 Q 
2"Mg 0-5 Small (<10) 


Our new data suggest strongly that for the nuclei near ?’Al itself, not 
even one heavy particle has been emitted. Thus, we would expect ?°Al 
formed by the simple ‘ stripping’ of a neutron from !4N, to be a very 
probable product which indeed it is. The yield of ?’Mg is barely one third 
of one per cent of 28Al and though loss of a neutron would doubtless be 
more probable from an excited state of ?°Al than loss of a proton, the 
difference of excitation for evaporation of these particles required is small 
and it is not easy to believe that so little ?’Mg would be found if much of 
the 28Al first formed were to break down. The higher yield of ?4Na 
(compared with 27Mg) shows even more conclusively that most of this has 
not been derived from 28Al as an intermediate stage, since an alpha 
particle requires more energy and faces a higher potential barrier than 
does a proton when leaving ?8Al. The production of *4Na thus requires 
some other process. We postulate one in which the “N abstracts 
nucleons, i.e. two protons and one neutron, from the “Al. Such a process 
would not appear inherently very likely, but it may be helpful to regard it 
as an interchange of a neutron and an alpha particle. Stated in this way 
it becomes a simple extension of the ‘ buckshot ’ principle. 

We shall now go on to consider what can be learnt from the shapes of 
the distribution curves of figs. 2 (a) and 2 (6). Since the energy distribution 
of the beam is continuous, as shown in fig. 1, these shapes do not represent 
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the excitation functions of the various products. They depend on these 
in part but are mainly determined by the other factors considered below. 


5.2. Distribution of the Heavier Products, from *4Na to 8K 


A striking feature of the distribution of these products is the rapid rise 
in yield, by a factor of the order 10, on passing from the cover foil to a 
maximum occurring at a depth usually between 10 and 20 from the 
surface. (In figs. 2 (a) and 2 (b) we have plotted the relative yields at 
depths corresponding to the middle of the foils, although this involves 
slight inaccuracy when the yields are changing rapidly with depth.) 
This initial rapid rise is clearly connected with the recoil of active products 
following the impact of the 14N ion. The depth in the target of the maxi- 
mum is of course determined in part by the recoil, but other factors appear 
to come into play and it is not easy to find convincing correlations between 
the depths of the maxima and, say, the mass number of the products. 

The next property of the curves to be discussed is their rate of fall-off 
in the deeper parts of the foil stack. In an earlier paper it was stated that the 
proportions of the different products varied by less than a factor of two at 
different depths. After the first 5, which was not then examined, this is 
indeed true down to at least 40. In the earlier work the 14N beam was 
not intense enough to give useful activities much below this depth. 
It is now clear, however, that there is a real change in the relative yields 
in the deeper foils. The final slope of the curves is unfortunately rather 
sensitive to the initial spectrum of 14N energies, which, as has been noted, 
varies somewhat from one bombardment to another. It is therefore 
unprofitable to attempt a quantitative discussion of minor differences 
observed. Some important general observations may however be made. 

Firstly, the curves show similarity in that they all approach asymptoti- 
cally to an exponential decrease. This is presumably connected with the 
fact that the range-spectrum of the incident 14N beam also shows a roughly 
exponential fall-off at the greatest ranges ; however, recoil phenomena will 
again come into play here. Secondly, the rate of fall-off is not always the 
same. We may note in this connection that, had the products all been 
derived by particle evaporation from the compound nucleus 4!Ca, we 
should reasonably expect a bias in favour of the lighter products in the 
upper foils, since there more reactions can occur in which energy is avail- 
able to cause evaporation of the many particles needed. But in fact the 
curves tend to show the opposite bias ; thus those for ?4Na, 27Mg, 28A1 
and *?P fall off less steeply than that for 34Cl. It is true that the fall-off 
for *°K is even less steep than the others, but in this case it is clear that 
8K must be formed by the fusion of the ?’Al with nearly all, if not quite all, 
of the N with subsequent particle evaporation. In either case it is 
therefore to be expected that an appreciable fraction of 38K nuclei formed 
in the early foils will still have sufficient excitation to evaporate an alpha 
particle. This process will diminish the *8K yield in this region and con- 
tribute to the *4Cl yield. The third observation we can make is that the 
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curves do not approach the exponential fall-off at equal rates. The 
remarks above give an immediate explanation of why the *4Cl curve 
should so quickly approach its asymptote (due to enhanced yields in say 
the first 40 1) while the 38K curve should do so slowly. 

The enhanced *4C] yields at the expense of ®8K is in effect a competition 
phenomenon rather analogous to the well-known case of excitation 
functions of processes such as («,n), («, 2n) ete. As a test of the purely 
“evaporation ’ model one would like to search for convincing cases of 
competition between processes leading to the lighter products. Such 
effects are however not at all obvious, although this might be because so 
many isotopes happen to be unobservable. Finally we notice that the 
28Al curve approaches its asymptote most rapidly of all. This fits very 
well with the supposition that ?°Al is formed simply by a neutron exchange, 
or stripping reaction, in which there is very little recoil in any direction. 


5.3. Distribution of the Light Products 


In our previous paper (Chackett and Fremlin 1954) we concluded that 
the N ions must be produced by a ‘ stripping’ process rather than by 
evaporation or fission. The evidence for this view was that the #=N atoms 
are observed in greatest yield at a considerable depth (about 60) in the 
foil stack, showing a very strong forward displacement from the site of 
formation. It is now clear from fig. 2 (b) that the mean forward range of 
the ®N atoms is even greater than that of the 14N ions from which they 
are derived. This can only mean that nearly all the product ®N nuclei 
must be projected in a very narrow cone. Since on kinematical grounds 
one would not expect such a strongly peaked angular distribution for 
incident 14N ions of low energy, it is probable that the cross section for the 
13N production rises fairly sharply with energy. In classical terms this is 
to say that the odd neutron in the 14N can most easily be removed, without 
disrupting the residual light nucleus, in a high-velocity glancing collision. 
The other nucleus so formed will of course be ?8Al, and this will have very 
little recoil. In this connection it should be remembered that *%Al is 
formed far more abundantly than is #N. It is therefore likely that the 
bulk of the 28Al corresponds to interactions of the 14N with the ?7Al, 
in which 18N as first formed is sufficiently excited to undergo break-down 
to lighter nuclei and so escape detection. It is nevertheless still true since 
only a neutron is transferred from one nucleus to the other, that the ?8Al 
suffers only a slight recoil in the direction of the beam. In the cases 
where the 14N and 27Al nuclei undergo still closer interactions, involving 
a larger transfer of momentum, the predominant reactions will be those 
leading to other products. | 

The distribution curve for 11C is noticeably broader than that for UN, 
This may be interpreted as indicating a larger interaction enetey in its 
production, which would lead to a wider angular distribution and nents to 
a greater variability of penetration for any given energy of the bombarding 


nucleus. Nevertheless the interaction energy must still be fairly small 
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even though three nucleons are removed from the 4N. The somewhat 
similar curve for 18F also implies a reaction mechanism involving only a 
moderate interaction energy ; here the 14N nucleus is evidently picking up 
nucleons. Some momentum will be carried by the remaining part of the 
27Al nucleus, but even if none were carried thus the velocity of the *F 
nucleus would be only 14/18 of that of the incident 14N, so that an earlier 
maximum in the distribution is naturally expected. Also it seems likely 
that the 18F ions have higher average state of ionization that the '4N, so 
that even for a given initial velocity their ranges would be less. Little can 
be said as regards the 10. This isotope is formed in too small a yield in 
the upper foils to enable us to distinguish it from other short lived products. 
It is however the principal product of those with half-lives of the order 
2-3 minutes at depths of 60-70. There can be no doubt therefore that 
150 also fits well into this group of products as regards its depth distri- 
bution and hence mode of formaticn. 

To conclude this section we should like to mention a few exploratory 
experiments on the distribution of products from 16O bombardments of 
aluminium. While there is some difficulty in correlating the results in 
detail which is probably connected with the lack of reproducibility of 
operating conditions of the cyclotron, the same general features are obser- 
ved. Thus 18N and 18F show.strong forward displacements indicative of 
formation by a stripping or a pick-up process ; the main heavier products 
(2°P, 34Cl, 38K) show distributions which are qualitatively similar to those 
found in the 14N bombardments. 


§6. CoNCLUSION 


The results clearly show that there are two main groups of products 
from 14N bombardments of ?’Al. These comprise a light group (14C, 4N, 
150 and !8F) which are formed in small yield by ‘ stripping ’ or ‘ pick-up ’ 
processes, and a heavy group (24Na, 27Mg, 28Al, °°P, 218i, 82P, 34C], 38K) 
mostly formed in much higher yield by the absorption by the ?*Al of part, 
or all, of the incident 14N ions. Similar results are obtained with 16O 
bombardments. 

Although they are not discussed in previous sections, we should like to 
add the following comments. 


(a) The isotopes of Ca and Se which might be produced by 14N and 
16Q bombardments of Al are not detectable by our techniques. 

() Many reactions will give rise to stable, as opposed to radioactive, 
products. 

(c) The lower relative yield of the ‘light group ’, %C, 18N, 150 and 18F 
may be due in part to the fact that, with quite moderate excitations they 
may split up into stable fragments. 

(d) The absolute yield of the products is not known accurately ; the 
total observable reaction cross section is a few per cent of the geometrical 
cross section. 
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(e) The reactions (14N, xn) and (260, wn) which are so well known with 
the heavier elements, might be invoked as an alternative mechanism for 
the production of say 38K, regarded as the daughter of a short lived 38Ca. 
But the value of x becomes impossible large for the production of even 
34C] from °4Ca let alone ®°P and 18F, and of course this mechanism cannot 
contribute to the production of any neutron-excess product. 
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Il. Inelastic Neutron Scattering in Manganese 


By Joan M. FREEMAN 
Atomic Energy Research Establishment, Harwell* 


[Received August 20, 1954] 


SUMMARY 


Gamma-rays of energies 126, 853 and 980 kev have been observed 
following inelastic neutron scattering in °°Mn, with neutron bombarding 
energies up to 1:3Mev. The evidence indicates levels in °°Mn at 
126-+3 kev and 979-L8 kev, the second level decaying mainly (~95%%) by 
cascade through the lower level. 


In order to study low excited states of °°Mn the gamma rays following 
inelastic neutron scattering in ®>>Mn have been examined, with both stilbene 
and sodium iodide crystal detectors, for neutron bombarding energies up 
to 1:3. mey. Manganese chips (1050 g) were enclosed in a thin cardboard 
container of annular form placed to receive neutrons emitted from the 
T (pn) reaction at a mean angle of 7-4° to the proton beam direction. The 
crystal detector was situated at the centre of the ring, on the axis defined 
by the proton beam, so that it recorded y-rays emitted from the manganese 
at a mean angle of about 83° to the neutron flux direction. The detector 
was shielded from direct neutrons by a polythene cone. Calibrated BF, 
counters were used to measure the neutron flux. 

At each neutron bombarding energy the pulse-height distribution from 
the crysta .was recorded on a thirty-channel kicksorter, first with the 
manganese ring in position and then with it removed. The two distri- 
butions were then subtracted to give the difference spectrum due to 
reactions in the manganese. Full details of the experimental method will 
be presented in another paper (Freeman, Lane and Rose 1955). 

The reaction was first examined with a 1 in. cube stilbene crystal and, 
for neutron energies greater than 1 Mev, Compton spectra corresponding 
to y-rays of energy about 850 kev were observed. Since iron is known to 
give y-rays of about this energy, following inelastic scattering (see, for 
example, Rose and Freeman 1953, Kiehn and Goodman 1954) and might 
possibly have caused the observed y-rays, the spectra due to an iron 
scatterer were measured at the same time. Examples of the difference 
spectra obtained for the two scatterers with 1-05 Mev neutrons are given 
in the inset of fig. 1, together with a calibration spectrum due to 47Cg 
(662 kev). A small difference in energy for the two y-rays was established. 
The Cs source and a ®Co source (1-17 and 1:33 Mev) were used for absolute 
y-ray energy calibration. The averages of a number of runs gave the 
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results: 848+15 kev for the y-rays from ®°Mn and 833+10 kev for the 
y-rays from **Fe. 

The relative y-ray yields at different neutron energies were obtained 
_ from the Compton spectra after corrections had been applied for the 
elastically scattered neutrons (which produce the steep rise below the 
Compton peak in fig. 1, inset), and for the backgrounds due to harder 
y-rays. The resulting excitation functions for Mn and Fe are shown in 
fig. 1. The curve for iron is similar to that obtained by Kiehn and 
Goodman (1954). The curve for manganese shows a threshold just below 
1 Mev, indicating a level in *>Mn at about 980-20 kev and showing that 
the 848 kev y-rays arise from a cascade process through a lower level. 
The absolute cross section values given in the figure were obtained by 
calculating the crystal efficiency for the detection of 840 kev y-rays from 
each scatterer, the calculation being based on absolute measurements 
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made previously with a standard 7Be (478 kev) source (see omen et al. 
1955). An assumption of isotropy of the y-rays (observed at ~90 ) Mie 
made. The standard deviation of the absolute cross section scale Is 40%. 
The investigation of manganese was continued with Byeocrunt iodide 
crystal (14in. diameter by 1 in.). For neutron energies ae e, the 
1 mev threshold, 850 kev y-rays were again observed (see further asap 
below) and 126 kev y-rays were also detected. The latter group was in 
present when the neutron energy was reduced below the threshold for the 
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850 kev y-rays, showing that in the cascade transition the 850 kev y-rays 
must precede the 126 kev y-rays. A typical 126 kev spectrum, obtained 
with 900 kev neutrons by subtracting the pulse-height distributions with 
and without the manganese ring in position, is given in fig. 2, together with 
calibration peaks due to americium and indium (y-rays 59-8 and 191 kev 
respectively, from Hollander e¢ al. (1953)). 

The y-ray measurements are summarized in the table, and indicate levels 
in 55Mn at 126-+3 kev and 979+8 kev. These results are in agreement 
with the observations of Hausman et al. (1952), who detected inelastic 
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Sodium iodide pulse-height distributions due to y-rays from Mn (nn’) at a 
neutron energy of 900 kev, and from Am and In. 


Observed gamma rays 


Nuclear 
levels 
Stilbene Nal Mean 
kev kev kev kev 
126+3 126525 126+3 
55Mn 848 +15 854-48 853 +7 979 +8 
~980 
56Fe 833-410 843-45 | 841.5 8414.5 


proton scattering to a 130 kev level and to a 1-0 mev level, and of Temmer 
and Heydenburg (1954) who observed 128 kev y-rays from manganese as a 
result of Coulomb excitation. 
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The y-ray spectrum from manganese for y-rays in the energy region 850 
to 1000 kev, with neutron energies from 1-1 to 1-3 kev, was studied with 
the sodium iodide crystal in a search for the ground state transition from 
the second (979 kev) level. For this purpose pulse-height distributions 
were obtained with manganese, iron and graphite scatterers. The latter 
was chosen because carbon is known to give no y-rays from inelastic 
scattering at the neutron energies used, so that the pulse height distri- 
bution with this scatterer represented the background due to y-emitting 
processes in the sodium iodide crystal, produced by scattered neutrons. 
The counting rate at the 850 kev peak in the manganese distribution was 
about twice the rate at the same pulse height with carbon, and for iron the 
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ratio was about five. At pulse heights a little above that corresponding 
to 1 Mev in energy the rates in all three distributions were approximately 
the same. The carbon distribution was fitted to the manganese distri- 
bution by normalizing to the same average number of counts in the ep 
channels (above 1 Mev) and was then subtracted from it to ane the 
difference spectrum due to manganese y-rays. The ys ee 
spectrum was obtained in the same way. The average eae te ce - 
sets of runs are given by the upper curves In fig. 3. The ee cor ae te i 
by reference to 1°’Cs and ©°Co sources, to y-ray energies of 854-L8 kev for 
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Mn and 843-5 kev for Fe. The manganese spectrum shows some sign 
of a small higher energy component. ‘To analyse this the Fe spectrum was 
normalized by the appropriate pulse-height and intensity factors so that its 
peak fitted the main Mn peak, and the two distributions were then 
subtracted. The result is shown by the lower curve in fig. 3. It can be 
seen that a small peak occurs at the position expected for 979 kev y-rays. 
The rate for coincident absorption of both 853 and 126 kev y-rays in the 
crystal was calculated from the crystal efficiency for 126 kev y-rays and 
was found to account for not more than 20% of the observed 979 kev 
peak. From this result and the relative peak intensities we conclude that 
the relative probability of emission of 853 and 979 kev y-rays from the 
979 kev level in Mn is about 20: 1. 

The work of Temmer and Heydenburg (1954) on Coulomb excitation of 
the first excited state (126 kev) of manganese demonstrates H, excitation 
and therefore the same parity for both the ground and 126 kev states. 
This result shows that the two competing y-rays from the second (979 kev) 
level must have the same parity ; the observed preference for the cascade 
transition then indicates that the order of multipolarity of the 979 kev 
y-ray is one unit higher than that of the 853 kev y-ray. 


I am grateful to Mr. B. Rose for valuable suggestions and criticisms 
during the course of this work. 


REFERENCES 


Freeman, J. M., Lanz, A. M., and Ross, B., 1955, Phil. Mag. (in the press). 

Hausman, H. J., Arten, A. J., Antuur, J. S., Benpsr, R. S., and McDotz, 
C. J., 1952, Phys: Rev., 88, 1296. 

Horianver, J. M., Pertman, I., and Sxazora, G. T., 1953, Rev. Mod. Phys., 
25, 469. 

Kinun, R. M., and Goopman, C., 1954, Phys. Rev., 93, 177. 

Rose, B., and Freeman, J. M., 1953, Proc. Phys. Soc. A, 66, 120. 

Temmer, G. M., and Heyprensure, N. P., 1954, Phys. Rev., 98, 351 and 
International Conference, Glasgow. 


bee liae 


Ill. Inelastic Neutron Scattering in Lithium Seven 


By Joan M. Freeman, A. M. Lanz and B. Ross 
Atomic Energy Research Establishment, Harwell} 


[Received August 25, 1954] 


ABSTRACT 


The cross section for the reaction 7Li(nn’)7Li*, in which the residual 
- nucleus is left in its first excited state (478 kev), has been measured as a 
function of neutron bombarding energy from the threshold up to 1-65 
Mev. The reaction was studied by observation, with both stilbene and 
sodium iodide crystal detectors, of the y-rays from the 7Li*. Details 
of the experimental method are given. The shape of the excitation 
curve for the reaction immediately above threshold indicates a strong 
S-wave contribution due to a (1—) level in the *Li compound nucleus. 
Above this region the curve rises to a peak at a neutron energy of about 
1:35 Mev, corresponding to a SLi level at 3-2 Mev (I’~1 Mev), which can 
be interpreted in terms of theoretical models as either a (1+) or, possibly, 
a (1—) state. 


§ 1. INTRODUCTION 


A REACTION of the type (nn’), in which the bombarded nucleus is left in 
an excited state after the emission of a neutron of energy lower than the 
bombarding energy, is generally studied by observation either of the 
inelastically scattered neutron or of the gamma-ray subsequently emitted 
as the residual nucleus reverts to the ground state (see, for example, 
Poole (1953), Kiehn and Goodman (1954), Eliot e¢ al. (1954)). 

Studies at a fixed neutron bombarding energy give data on the energy 
levels of the residual nucleus, and cross sections. Additional information 
can be obtained by varying the neutron energy and measuring the 
threshold and excitation function for a particular reaction. 

The gamma-rays following inelastic scattering are usually observed 
with a sodium iodide crystal but can also be detected with an organic 
scintillator (Rose and Freeman (1953)). We have used both types of 
crystal for the observation of the gamma-rays and have found that the 
organic phosphor has particular advantages in the measurement of the 
excitation function and cross section in the neighbourhood of the 
threshold. A sodium iodide crystal is more useful for higher neutron 

arding energies. 
pe eee, ae use of the two types of detector will be given in the 
next section and the experimental details of the study of the reaction 


7Li (n n’) 7Li* will then be described. 


+ Communicated by the Authors. 
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§ 2. Tue CHOICE OF THE SCINTILLATOR FOR y-RAY DETECTION 


A crystal detector set up to record y-radiation following inelastic 
neutron scattering will also receive neutrons scattered elastically and 
inelastically from the target material. It is therefore necessary to be 
able to distinguish between pulses in the detector due to the target 
y-rays and pulses due to the products of neutron reactions in the crystal. 
In a sodium iodide crystal good resolution and efficiency are obtained 
for the y-ray lines, but the neutrons induce inelastic scattering and 
capture reactions in the sodium and iodine nuclei, with consequent 
y-ray emission. The latter constitutes an appreciable and varying 
background, so that near the threshold for the reaction under investiga- 
tion, or whenever the cross section is small, the y-ray peak due to this 
reaction may be less intense than the background groups. In this case 
quantitative estimates and even identifications are difficult. For 
example, when the reaction “*Li(nn’) 7Li* was studied with a Nal 
crystal, the 478 kev y-rays from 7Li* were found to be superimposed 
on a broad group at about 430 kev, due to inelastic scattering in sodium 
and iodine (see fig. 7, § 3.3). On this account quantitative yield measure- 
ments could not be made near the threshold. 

In an organic crystal the y-rays are detected from the Compton 
recoils and the main effect of scattered neutrons is the production of 
proton recoils, although some capture also occurs. The response of 
organic phosphors to protons is non-linear (Birks 1951), and in the energy 
region of interest in the study of the lowest excited states of nuclei, the 
proton recoil pulses due to neutrons are much smaller than the electron 
recoil pulses due to y-rays of about the same energy. The y-ray groups 
are then easily distinguishable from the scattered neutron groups. As 
far as the y-ray background due to neutron capture is concerned, the 
Compton efficiency of the crystal is greater in the energy range of interest 
to us than for the unwanted harder y-rays and this background is found 
to be small enough to require only a small correction. 

The resolution obtainable with y-ray spectra in an organic crystal, 
although not as good as that of a good Nal crystal, is better than is 
perhaps generally realized. Some typical spectra obtained with a one 
inch cube stilbene crystal and an E.M.I. Type 6260 photomultiplier are 
shown in fig. 1. The pulse-height distributions were recorded on a 
thirty-channel kicksorter. Curve A gives the spectrum due to the 
1-277 Mev y-rays from **Na; the spectrum shows the characteristic 
peak, and a well defined end point is obtained by extrapolation of the 
linear part of the cut-off. Curve B is the distribution due to the 1-171 
Mev and 1-331 Mev y-rays from ®°Co. The peaks due to the two y-rays, 
which are 13%, apart in energy, are almost resolved. The end point 
corresponds, by reference to the 7*Na and other calibration sources, to 
1:33 Mev, and if the *Na spectrum shape is used to subtract away the 
higher energy component, the residual spectrum (curve C) corresponds 
to the 1-17 Mev component. Curve D is the spectrum due to the 
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y-radiation from the decay of **Mn, which gives y-rays of 840 kev, 1-81 
Mev and 2-06 Mev in the intensity ratio 10:3:2. The presence of two 
lines in the higher energy part of the spectrum is apparent from the 
change of slope, and the energies of the two components can be fixed 
by comparison with other calibration spectra. 

Recoil proton spectra obtained from monoenergetic neutrons in the 
same stilbene crystal are shown in fig. 2. Examples for three neutron 
energies are given ; in each case the energy spread of the neutrons was 
~20 kev. The broken curve shows, on the same pulse-height scale, 
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the spectrum due to the 662 kev y-rays from 1°7Cs, and illustrates the 
difference in response of the crystal to neutrons and y-rays, as mentioned 
above. peer ne 
The stilbene crystal calibration for both y-rays and. neutt ons is given in 
fig. 3, the extrapolated end points of the spectra being plotted against 
sii e trons. 
ry energy of the y-rays and neu . 
P The Gees crystel is thus very suitable for measuring the ys adiation 
resulting from inelastic neutron scattering for neutron energies 1 ate 
neighbourhood of the threshold for excitation of a given level. The 
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method can be extended up to neutron energies for which the spectrum 
produced by the elastically scattered neutrons begins to overtake the 
Compton peak of the y-ray group. If the y-ray transition takes place 
directly from the excited level under investigation to the ground state, 
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it is possible for the reaction to be followed for neutron energies up to 
somewhat more than twice the threshold energy. If the y-radiation 
oceurs by cascade through a lower level, the observable range is corre- 
_spondingly reduced. We have found, as will be described below, that 
there is adequate overlap between the neutron energy ranges for which 
the organic and Nal crystal detectors give good quantitative results. 


§ 3. THE Reaction 7Li (nn’) 7Li* 
3.1. Hxperimental Method 


The arrangement used for the study of the reaction ? Li (n n’) 7Li* is 
shown schematically in fig. 4. The source of neutrons was the T (pn) 
3He reaction produced by protons from the Harwell Van de Graaff 
machine bombarding a zirconium tritide target (~ 200 ug/cm? of Zr). 


Fig. 4 
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Experimental arrangement (schematic). 


The lithium metal scatterer, which had been cast in a steel container 
and then removed, dipped in vaseline and wrapped with cellotape, was 
in the form of an annular ring of outside diameter 11-5 cm, inside diameter 
5-6 cm and thickness 4:2 cm. It was placed with its axis in the direction 
of the proton beam and its centre 33 cm from the target spot so that it 
received neutrons emitted at a mean angle of 7-4° to the proton beam 
direction. The 1 in. cube stilbene crystal detector was mounted centrally 
inside the lithium ring, and a polythene cone 18 em long set up 
between the source and the crystal served to shield the latter from 
neutrons emitted directly from the ZrT target. The energy of neutrons 
bombarding the lithium ring was calculated by reference to the threshold 
for the T (pn) reaction (1-02 mev). The energy spread of the neutrons 


reaching the ring was 20 kev. 
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A boron trifluoride long counter embedded in paraffin wax was placed 
to count neutrons emitted at a mean angle of 60° to the proton beam 
direction and was used as a monitor throughout the runs ; it was after- 
wards compared with a second absolutely calibrated BF; chamber set 
up at a mean angle of 7:-4° to the proton beam, so that the neutron flux 
in the lithium ring could be determined. 

At each neutron bombarding energy the spectrum of pulses produced 
in the crystal detector was recorded first with the lithium scatterer in 
position and secondly with it removed. The*difference between the two 
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Stilbene crystal spectra with 720 kev neutrons bombarding lithium. Curve A: | 
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spectra, normalized to the same number of neutron monitor counts, then 
represented the pulse distribution due to processes in the lithium ring. 
The results of a typical pair of runs for a neutron energy of 720 kev are 
given as log-log plots in fig. 5. Curves A and B are the pulse distributions 
obtained with and without the lithium ring respectively and the broken 
curve C, with points marked as circles, is the difference spectrum. The 
steeply falling lower energy part of curve C is due mainly to elastically 
scattered neutrons, and the group at higher energy, showing a Compton 
peak at a pulse height of about 18-5 volts, is due to 478 kev y-radiation 
from the first excited state of “Li formed in the (nn’) reaction. For 
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comparison the spectrum obtained with a 7Be source, which also produces 
y-rays from the first excited state of 7Li, is shown by the full curve D in 
fig. 5. 
; 3.2. The Excitation Function near Threshold 

Some of the difference spectra, obtained at the neutron bombarding 
energies indicated, are shown in a linear plot in fig. 6. The pulse distri- 
bution due to 478 kev y-radiation, with an end point at 26-5 volts, is 
evident in the curves taken above the threshold (#,—546 kev) for the 
inelastic process. All the curves, including that taken below the thres- 
hold, indicate the presence of a small amount of harder y-radiation. 
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Spectra due to y-rays from “Li (n n’) at the neutron energies marked. 


This was found, from measurements of the end point, to have an energy 
of about 2-2 mev, and was assumed to be due mainly to capture by 
protons in the stilbene crystal of neutrons scattered from the lithium 
ring; a small contribution could have been due to slow neutrons fon 
the polythene scatterer being captured by “Li. The shape of a spec- 
trum due to this harder y-ray background was assumed to remain the 
same over the range of neutron energies used and to be represenved by 
the average of all the difference spectra taken below es . n threshold. 
The magnitude of the correction to be applied to each of the distri Bitigns 
obtained above threshold was then determined by fitting this average 
in each case at points above the 478 kev end point. es corrections 
were found to be fairly constant for all neutron bombarding energies. 
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Before estimations of the relative 478 kev y-ray yields could be made, 
the effect had also to be calculated of the contribution to the spectrum 
due to the high energy tail of the elastically scattered neutron group. 
The latter was obtained for each neutron bombarding energy #,, as 
follows. 

From the geometry of the system, the scattered neutrons which could 
be received by the crystal had a mean energy of 3/4 H,, (corresponding to 
90° scattering) and a spread of +0-:07 #,,. The shape of the spectrum 
expected for this group was calculated from the appropriate distribution 
obtained with monokinetic neutrons (fig. 2). The curves derived in this 
way were found to fit the difference spectra taken below the (nn’) 
threshold. The use of log—log plots (as in fig. 5), which could be super- 
imposed, was found useful for this curve fitting. The derived curves 
were then fitted to the experimental curves for neutron energies above 
the threshold to determine in each case the pulse height at which the 
scattered neutron contribution to the 478 kev y-spectrum became 
negligible. The counts in each spectrum were integrated from this 
pulse height up to 23 volts. Use was then made of the known shape of 
the 478 kev distribution, obtained from a 7Be source, to normalize these 
integrations to the same pulse height interval. From the values of the 
integrated gamma-ray counts obtained at different neutron energies and 
the corresponding neutron flux measurements, the relative excitation 
function for the reaction was thus obtained. The results are given by 
the points marked as circles in fig. 8, the standard deviations being 
indicated for each point. Most of the points represent the mean of two 
or three separate measurements. The greater errors at the higher neutron 
energies are due to the fact that the scattered neutron groups have reduced 
the usable ranges of the gamma-ray spectra. 

A small correction to the excitation function had to be made for the 
loss of neutron flux in the lithium ring due to scattering. The magnitude 
of the total scattering cross section was obtained from published data 
(Bockelman et al. 1951). For neutrons of initial energy less than about 
700 kev the scattered neutrons were effectively lost, but for higher initial 
energies some of the elastically scattered neutrons were still energetic 
enough to be able to contribute (at a reduced rate because of their loss 
of energy) to the inelastic scattering reaction. Appropriate corrections, 
which amounted to a few per cent, were calculated for this effect. 

The higher energy points and absolute cross section values given in 
fig. 8 will be discussed in the next two sections. 


3.3. The Excitation Function at Higher Neutron Energies 


An attempt was made to record the inelastically scattered neutrons 
at higher energies, but the energy spread of these groups was too great, 
due to the geometry of the system, to allow satisfactory measurements. 
The stilbene crystal was then replaced by a 14 in. diameter x1 in. thick 
Nal crystal for examining the “Li (nn’) excitation function above 1-1 
Mev. An example of a difference spectrum obtained by subtracting runs 
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with and without the lithium scatterer in position is shown by curve A 
in fig. 7 ; the neutron energy was 1:22 Mev. The corresponding difference 
spectrum obtained with a graphite scatterer of the same dimensions is 
given by curve B. Since carbon cannot give y-rays from inelastic scatter- 
ing at this neutron bombarding energy, the spectrum B must be due to 
_ neutron-induced reactions in the Nal crystal. The shape of this spectrum 
was found to be insensitive to small changes of neutron energy, so that 
it could be assumed to represent the shape of the background component 
in the spectrum obtained with the lithium scatterer. The carbon curve 
B was therefore fitted to the lithium curve A at points well above the 
478 kev peak and the difference was taken to give the spectrum due to 
the 478 kev y-rays from the 7Li (nn’) reaction. The curve C in fig. 7 
shows this difference which agrees well in shape and position with cali- 
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NaI spectra due to lithium (curve A) and carbon (curve B) for 1-22 Mev 
neutrons. Curve C is the difference spectrum. 


bration spectra obtained with a 7Be source. The relative y-ray yields at 
different neutron energies were obtained by integrating over the peaks 
of curves of type C. The neutron flux estimates were again corrected 
for the effect of neutrons elastically scattered in the lithium ring. HIN 
yield results are given by the points marked as crosses in fig. 8. It was 
found possible to extend the Nal crystal points with peancuaule oy 
down to a neutron energy of 800 kev, thus allowing enough overlay ae 
fitting to the stilbene crystal results. This was done by deawing 4 io 
ve through the stilbene points, finding the ratios of Nal points to 
curve and then taking the meai ratio tor 

The standard deviation of the mean ratio 


best cur 
corresponding points on this 
adjusting all the Nal results. 


was 5%. 
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3.4. The Absolute Cross Section ; 

Two partly independent determinations of the absolute cross section 
for the 7Li (n n’) reaction were made, one for results obtained with the 
stilbene crystal at a neutron energy of 712 kev and the other for Nal 
crystal measurements at 1320 kev. 

The efficiencies of both the stilbene and Nal crystals for the detection 
of 478 kev y-rays were determined with the aid of a small 7Be source of 
known strength. In each case the total counting rates in a fixed pulse 
height interval of the y-ray spectrum were measured as the source was 
moved to different points in the space which the lithium ring occupied 
during the scattering experiment. The pulse height interval was chosen 
to be the same as that used in the integration of the counts due to y-rays 
from the (nn’) reaction as described in §§ 3.2 and 3.3. The absorption 
due to the total lithium ring thickness was then measured and found 
to be about 10% in agreement with the theoretical value. It was thus 
possible to calculate the counting rate expected with each crystal for a 
478 kev y-ray source of known strength distributed uniformly throughout 
the lithium ring. The y-rays from the first excited state of lithium were 
assumed to be isotropic (Ajzenberg and Lauritsen 1952). 

The absolute neutron fluxes in the scatterer were obtained from the 
calibrated neutron monitors and the known geometry of the system, 
corrections being applied for the loss due to elastic scattering as discussed 
in § 3.2. The number of *Li atoms per cm? in the ring was estimated 
from its mass and volume after removal from the steel case in which it 
was cast. Allowance was made for the isotopic concentration of ®Li 
CTD oe): 

The value of the cross section obtained with the stilbene crystal and 
712 kev neutrons was 0-050 10-74 cm?. The standard deviation for 
each factor in the cross section calculation was as follows: integrated 
counts from the pulse height distribution, 4% ; absolute crystal efficiency, 
12% ; effective neutron flux in the lithium ring, 13°; number of 7Li 
atoms per cm?, 10%. The total standard deviation was 21%. 

The cross section found with the NaI crystal and 1320 kev neutrons 
was 0:23 x 10-*4 cm? with the same standard error. From the excitation 
curve of fig. 8 the ratio of the yields at 1320 kev and 712 kev is 3-81 
(+5%), so that the stilbene measurement of 0-050 x 10-24 em? at 712 kev 
is equivalent to a cross section of 0:19 10-24 cm? at 1320 kev. We take 
therefore the average value of the cross section for the point at 1320 kev 
to be 0-21 x 10-24 cm? with a standard deviation of 19°%%. The ordinates 
given in fig. 8 are based on this value. 


§ 4. Discussion 


We now discuss the analysis of the excitation curve in fig. 8. The 
trend of the experimental points near the bombarding energy 1:35 Mey 
enables us to assume the presence of a resonance level in ®Li centred at 
this energy and with a total width of roughly 1 Mev. The existence of 
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this level can also be inferred from the measurements of the total cross 
section of “Li (Bockelman ef al. 1951). 

The first problem in fitting an excitation curve is to determine the 
predominant angular momenta involved in the reaction. Assuming that 
the ground and first excited states of 7Li have the same parity, we can 
immediately deduce that ingoing partial waves of one parity give rise to 
outgoing waves of the same parity. To try to discover which partial 
waves contribute most to the reaction, we first examine the energy 
dependence of the excitation curve near threshold. In fig. 8 (inset) 
curves A and B give the predicted dependences assuming s-waves in and 
out, and p-waves in and out, respectively (Wigner 1948). Higher 
partial waves give much steeper curves than B. It can be seen that curve 
A fits the observed points much better than curve B, not only at the 
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fit on the assumption (i) that the resonance is (1—). Curve D gives 
the best fit on the assumption (ii) that the resonance 1s (1+), super- 
imposed on a strong s-wave background given by curve A. An inter- 
action radius R of 4-20 10-" cm is assumed throughout. 


“threshold, but up to 120 kev above it. This implies the presence of a 
strong s-wave contribution to the reaction. Such a contribution must 
1—) levels of the compound nucleus (assuming that the two 


ise from . } 
aris ( (3/2—) and (1/2—) (Ajzenberg and Lauritsen 


states of 7Li have spins 


28 Joan M. Freeman, A. M. Lane and B. Rose on 


/ 


1952)). Ifwe assume that only one level at some energy #, and of reduced. 
widths yz and y,2, is responsible, our threshold fit and the measured 
absolute value of the cross section give the relation 

Yin * Yout 

a ee 

(# ,—-# Sip 
where a standard Breit-Wigner formula has been used and #, is the 
threshold energy (478 kev in centre-of-mass). It is interesting to note 
that an indication of strong s-wave levels in *Li has also been found in the 
analysis of the thermal neutron elastic scattering cross sections of “Li 
(Ajzenberg and Lauritsen 1952). 

It is tempting to assume that the level in our excitation curve at 1-35 
Mev is indeed the (1—) level responsible for these strong s-wave effects. 
However, before making such an assumption, one must consider if any 
other spin assignments to this level are plausible. We can immediately 
rule out any odd-parity assignments other than (1—) on the grounds that 
this would imply formation or decay of the compound state by d- or 
higher waves, which would lead to reduced widths far in excess of the 
Wigner limit (h?/MR?). An even parity state, excited by p-waves, is 
however possible, and the excitation curve could be interpreted as being 
due to a p-wave resonance at 1:35 Mev superimposed on an s-wave 
background (represented by curve A in fig. 8) due to a (1—) level at 
a higher neutron energy. This p-wave state could be either (1-+-) or (2+). 
Spin 0 is excluded because the observed cross section, after subtraction 
of the s-wave contribution, is greater than 7A?/8, the maximum allowed 
for this spin. We are therefore left with the following possibilities for 
the resonance level : 

(i) a (l—) state, with formation and decay mainly by s-waves, but 
with the possibility of some d-wave contributions, 

(ii) a (1+) or (2+-) state, with formation and decay by p-waves, the 
level being superimposed on the strong s-wave background given by 
curve A. 

This ambiguity could presumably be resolved by angular distribution 
measurements. In the absence of such data, we can take advantage of the 
recent advances in our theoretical ideas on nuclear structure. It has been 
established (Inglis 1953, Lane 1953, 1955) that an intermediate coupling 
shell-model gives an extraordinarily accurate description of the structure 
of light nuclei. The success of this model for states of equivalent nucleons 
(Inglis 1953) and its extension to mixed configurations (Lane 1954) often 
enables us to identify an energy level merely from its position in the 
spectrum. It is an unfortunate coincidence that such a model predicts 
both a (1—) and a (1+-) level in *Liin our range ofinvestigation. Further- 
more, we shall see below that detailed fittings of the experimental curve 
in the two cases give two sets of nuclear constants, each of which is 
consistent with the theoretical predictions. It seems however to be 
worth giving the details of the two cases to show how experimental 


material can be exploited in the light of recent advances in nuclear 
theory. 
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(i) Interpretation of the Resonance as a (1—) State mainly of the T'ype (1p)?2s 
(Lane 1954) 


An initial attempt to fit the whole excitation curve with only s-waves 
contributing resulted in poor agreement; either the theoretical Breit— 
Wigner curve fitted the points below 1-3 mev and fell below those above 
this energy, or vice-versa. The best fit (curve C in fig. 8) was obtained 
by introducing d-waves according to the parameters 

Eye3=1-35 Mev (lab.) 

Las (£.5)= 1-14 Mev (lab.) 

Yin(8)= Yout(8)= yin(d) = your(d) =12% #?/M R?. 
(The Thomas correction (Thomas 1951) was not taken into account here.) 
The fit was fairly insensitive to variations in the ratios of the reduced 
widths by factors up to 2, except for the requirement that yj(d) had to 
be at least as great as the mean of the two s-wave reduced widths. 
According to an analysis by one of us (Lane 1954), the ten or so levels 
in the Ip shell that are known to be formed by s-waves are all character- 
ized by having reduced widths near the single particle value 


h2 
(~ 0:8 wri) . 


This result implies that an incident s-nucleon couples so weakly 
to a target nucleus in the Ip shell that the latter is not appreciably 
excited above its ground state. However, in our case, the first excited. 
state of 7Li is so close to the ground state that even a quite weak coupling 
might be expected to excite it. This expectation is borne out by the 
near equality of the experimental values of the two s-wave reduced 
widths. The large d-wave reduced widths are rather surprising however, 
for they indicate a large amount of configuration mixing between the 
shell-model states (1p)? 2s and (1p)? 1d. There is no obvious sign of 
such a mixing in the other known levels formed by s-waves. It should 
be noted that the presence of these d-wave effects has to be assumed in 
our case not only to fit the excitation curve, but also to bring the sum of 
the four reduced widths of the compound level up to a value near the 
single-particle value and so to satisfy the sum rule (Wigner 1949). The 
two s-wave reduced widths alone can only account for about 30% of this 
figure, whereas the sum over all four reduced widths can amount to 70% 
(the Thomas correction tends to increase this latter value). 


(ii) Interpretation of the Resonance as the (1--) ('=1) State from the 
Configuration (1p)* (Ings 1953) 

Curve D in fig. 8 shows the result of fitting a Breit-Wigner curve to the 
experimental points on the assumption that the resonance is (1 +) (p- 
waves in and out), superimposed on the strong s-wave background given 
by curve A. Almost identical fits could be obtained. for aul ratios of the 
two reduced widths z=y2/yout from 1:10 to 10:1, provided that the 


resonance energy was suitably adjusted in the range 1-47 to 1:30 Mev. 
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A more severe restriction on the values of x is provided by the absolute 


p-wave cross section (0:14--0-04 barns at 1-3 Mev), which requires x to lie 
in the ranges 0-1 to 0-4 or 1:0 to 3-0. The fitting to the excitation curve 
demands a total width at 1:3 Mev of about 1-15 Mev, and this implies 
(for all x) that the sum of the two reduced widths is 0-7 (4?// R?) to within 
20%. 

Using the formulae given by Lane (1953) for reduced widths in terms of 
shell-model wave functions, we find that the predicted values of y;, and 
You, are 0-5 and 0:5 for the 38P[31] state in L-S coupling or 0-2 and 0-75 
for the state of J=1, T=1 in j-j coupling, the unit being the single 
particle lp reduced width (~ 0-7 (h?/MR?)). Both of these sets of values 
are consistent with the rough ranges of experimental values given above. 

We conclude from the above analysis of our data that the resonance 
level at 1-35 Mev could be either a (1—) or a (1+) state. Although both 
assignments are consistent with the data, the (1+-) assignment is slightly 
preferred because (a) a rather better fit is obtained with the excitation 
curve (compare curves C and D), and (b) the fitting of the (1—) level 
implies an unexpectedly large amount of configuration mixing between 
(1p)? 2s and (1p)? 1d. 

ACKNOWLEDGMENTS 


We should like to thank Dr. W. D. Allen for valuable discussions and 
interest in the work and for providing the calibrated BF, counters, 
Mr. W. B. Woollen for casting the lithium ring, Dr. R. C. Hanna for 
providing a calibrated “Be source, and Mr. A. E. Pyrah for his efficient 
running of the Van de Graaff machine. 


REFERENCES 


AJZENBERG, F., and Lauritsen, T., 1952, Rev. Mod. Phys., 24, 321. 

Brgs, J. B., 1951, Proc. Phys. Soc., 64, 874. 

Bocketman, C. K., Minuer, D. W., Apatr, R. K., and Barscuaty, H. H., 
1951, Phys. Rev., 84, 69. 

Exiot, E. A., Hicks, D., Becuran, L. E., and Hatpan, H., 1954, Phys. Rev., 
94, 144. 

Ineis, D. R., 1953, Rev. Mod. Phys., 25, 390. 

Krenn, R. M., and Goopman, C., 1954, Phys. Rev., 98, 177. 

Lane, A. M., 1953, Proc. Phys. Soc. A, 66, 977; 1954, A.E.R.E. (Harwell) 
Report T/R. 1289; 1955, Proc. Phys. Soc. (to be published) 

Pootg, M. J., 1953, Phil. Mag., 44, 1398. 

Ross, B., and Freeman, J. M., 1953, Proc. Phys. Soc. A, 66, 120. 

Tuomas, R. G., 1951, Phys. Rev., 81,, 148. 

Wiener, E. P., 1948, Phys. Rev., 78, 1001; 1949, Amer. Jour. Phys., 17, 99. 


Paelte 
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SUMMARY 

The supersonic flow past flat-nosed circular cylinders with annular 
gaps of rectangular section was investigated by schlieren photography and 
pressure measurements. The length of the gap was varied in small steps, 
from low values, where the gap formed a single dead-air region, to high 
values, where the flow in the gap consisted of two dead-air regions with 
an intermediate region of attached flow. In particular, the length of the 
gap for which the flow changed from the type with one to the type with 
two dead-air regions was determined. The experiments throw some 
light on the details of the flow in the dead-air regions, and it was found 
that the velocities inside these regions may be of considerable magnitude. 
Qualitative information was obtained about the effects of the length of 
the gap on the drag caused by the gap. 


§ 1. INTRODUCTION 
ConsipER a long circular cylinder of diameter D placed with its axis in the 
direction of a uniform supersonic stream. At a station far enough away 
from the front of the body for the static pressure to have reached its 
free-stream value, an annular ring of rectangular section is cut away from 
the cylinder, reducing its diameter to d over a length a. 
Fig. 1 


= 
as 


Flow with short gap. 


If the gap thus formed is short, as in fig. 1, the boundary layer which 
separates at the upstream edge of the gap would be expected to form 
a narrow almost straight mixing region passing across the gap in a direction 
nearly parallel to the axis of the cylinder. The gap itself will form a 
amas 
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so-called ‘dead-air region’ in which the velocities are small compared 
with the free-stream velocity, and the pressure is nearly equal to the free- 
stream static pressure. Air will be supplied to the dead-air region at the 
downstream end of the gap, where the mixing region meets the sharp 
corner of the body, and will be drawn away from the dead-air region along 
the mixing zone. 

If on the other hand the gap is long compared with the diameter of the 
cylinder, as in fig. 2, the boundary layer will no longer be separated right 
across the gap, but the flow pattern will be made up of three parts, namely 
two separated regions at the ends of the gap with an intermediate region of 
attached flow. 

At the upstream end of the gap the flow will be similar to that at the 
base of a body of revolution with a supporting sting. At the corner there 
is an expansion fan which bends the mixing region towards the axis, and at 


Fig. 2 


Flow with long gap. 


the bottom of the gap the mixing region is bent parallel to the axis by a 
shock wave or compression fan. The approximately conical dead-air 
region receives air at its downstream corner, and air is drawn out along 
the mixing region. 

This part of the flow which has been discussed by many authors (see 
for example Chapman 1950) will be followed by a region where the mixing 
zone is again attached to the body as a boundary layer. 

This boundary layer will separate under the influence of the down- 
stream end of the gap, and the flow here will be similar to that investigated 
by Mair (1952) who photographed the boundary-layer separation on probes 
in front of blunt-nosed bodies of revolution in supersonic flow. The 
dead-air region is again approximately conical, but the shock wave is 
now at the beginning of the mixing region, and the expansion at its end 
turns the flow in the direction of the axis of the body. The mechanism of 
the supply of air to and extraction of air from the dead-air region is 
similar to that at the upstream end of the gap. 

The experiments described in this paper were made to investigate the 
changes in flow pattern when the length @ of the gap was varied from a 
very small value where the flow pattern was of the first type described 
above, to a very large value where it was of the second type ; in particular, 
it was the aim to find out whether the transition from one type of flow 
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to the other was gradual or sudden. At first sight the photographs might 
seem to indicate that the change-over took place suddenly, but a more 
careful inspection combined with measurements of the pressure distri- 
bution along the gap showed that the sudden change was preceded by 
gradual changes as a approached its critical value. It was evident from 
the pressure measurements that the velocities in the so-called dead-air 
regions may have been of considerable magnitude. 

The main parameters of the problem are the geometry of the gap, i.e. 
D, d, and a, the free-stream Mach number M , and some Reynolds 
number based on the free-stream velocity and kinematic viscosity and on 
some typical length. This may be either the length L of the cylinder 
upstream of the gap or, perhaps better, the displacement thickness of 
the boundary layer immediately upstream of the gap. 

In the experiments, D, L, and the free-stream conditions were kept 
constant throughout, so that both the Reynolds numbers defined above 
were constant. Two values of d were used, and a was varied over a wide 
range. 

The experiments throw some light on the flow inside the dead-air 
regions and they also show that the drag caused by such gaps on bodies of 
revolution may be considerable even for small values of the length of the 
gap. 

§ 2. NoTaTION 
diameter of circular cylinder. 
length of cylinder upstream of gap. 
free-stream Mach number. 
pressure coefficient. 
free-stream velocity. 
length of gap. 
inside diameter of gap. 
static pressure. 
free-stream static pressure. 
distance along gap measured from its upstream end. 
free-stream density. 
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§ 3. EXPERIMENTS 


The experiments were made in an intermittent pce gree 
supersonic tunnel with a 4 in. x 5 10. working section. The stagnation 
conditions were approximately atmospheric, the free-stream Mach 
number M, was 1-965, and the Reynolds number was 34x 10° per inch. 
The flow was observed with a conventional two-murror schlieren system 
using a continuous mercury-vapour lamp for visual peer amane woe 
the photographs were taken either with a flash tube of effective dui ation 
about six microseconds or with a spark of duration just under one 
microsecond. The pressures were measured with absolute mercury 
We, ccodels are shown in fig. 3. They consisted of silver-steel rods 
of diameters d=} in. or } in. with a fixed and a moveable circular cylinder 
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both of outside diameter D=} in. and length L=1lin. The first cylinder 
was always fixed with its upstream end flush with the end of the rod, 
whereas the second cylinder could be fixed in any position, so that the 
length a of the gap could be varied continuously from 0 to 2inches. The 
steel rods were supported by cross-rods downstream of the working 
section. 

The length of the first cylinder was twice its diameter. According to 
experiments by Holder and Chinneck (1954) on flat-nosed cylinders at a 
Mach number of 1-82, the free-stream static pressure was reached at a 
distance of about one and a half diameters downstream of the flat nose. 
It can therefore be assumed that the pressure just upstream of the gap was 
equal to the free-stream static pressure in all the experiments described 
here. 


Fig. 3 


The models. 


For the pressure measurements the + in.-diameter rod was replaced by a 
1 in.-outside-diameter steel tube with twenty-one 0-015 in. holes drilled 
at equal intervals of 0-1 in. along the tube, the centre of the hole farthest 
upstream being in the plane of the base of the first cylinder. The holes 
were filled with soft solder and redrilled one at a time for the pressure 
measurements. 

The flows past both models were photographed for a number of values 
of a from 0-1 in. to 2-0in. It was found that to get as much information 
as possible, photographs with both flash-tube and spark were needed for 
each value ofa. A flash photograph gives a clear picture of the mean flow 
whereas a spark photograph shows details of the turbulence in the flow. 

Some typical photographs are reproduced in figs. 4-15 (Plates 1 and 2). 
They were all taken with the schlieren knife edge parallel to the direction 
of the free stream and arranged so that a positive density gradient upwards 
gives an increase in light on the photograph. Figures 4 and 9 (Plate 1) 
were taken with the flash-tube and the others with the spark discharge. 

The pressures were measured on the model with d=+ in. for values of a 
from 0-1 in. to 2-0 in. with a increasing in steps of 0-lin. Some typical 
pressure curves are shown in fig. 16 where each curve gives the pressure 
distribution across the gap for a particular value of a. The pressures are 
expressed as the non-dimensional pressure coefficient 

pee Seay ah 
2 PiU 
where p is the measured pressure, and p,, p;, and U, are the static pressure, 
the density, and the velocity in the free stream. 
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Fig. 16 


a= 0.6” 


Pressure distributions. 
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§ 4. DiscUSSION OF THE EXPERIMENTS 


From the visual and photographic observation of the flow pattern the 
value of a for which the flow changed character could be determined. 
For values of a below this critical value the flow was of the first type with 
one dead-air region extending across the whole of the gap. For values 
of a above the critical value there were two dead-air regions with an 
intermediate region of attached flow. 

For the body with d=1in. the critical value of a was found to be 
1:250in. This value of a was well defined. The flow was observed 
during a great number of tunnel runs for a=1-225in., 1-250 in., and 
1:275in. For a=1-225in. the flow was always of the type with one 
dead-air region, whilst for a=1-275 in. there were always two dead-air 
regions. For a=1-250 in. both types were possible, and either of them 
might persist throughout one tunnel run (about 30 sec), but in some runs 
a change-over was observed, and it was in all cases from flow with one to 
flow with two dead-air regions. 

For the body with d=tin. the critical value of a was found to be 
1-650 in., but this value was less well defined than for the body with 
d=jin. For arange of values of a around 1-650 in. the flow was unsteady, 
and photographs taken at random using the one-microsecond spark 
showed all intermediate stages between flows with one and two dead-air 
regions. Figures 11-14 (Plate 2) show photographs taken for a=1-650 in. 

It is believed that this difference in behaviour of the two bodies is not 
aerodynamically significant. The unsteadiness of the flow over the 
thinner body was probably caused by flutter of the model and might not 
have been present if it had been perfectly rigid. 

It is easy to see that the critical value of a must increase as d is decreased. 
The conical dead-air regions in the gap on one body have very nearly the 
same cone angles as the corresponding regions on another body with a 
different value of d, and the smaller d is the larger a therefore has to be to 
accommodate both regions, provided D is kept constant. 

Apart from the difference in the critical values of a, and the unsteadiness 
found on the body with d=} in., the flow patterns around the two bodies 
were very similar, and the detailed discussion will therefore be based on 
the body with d=jin. for which both photographs and_ pressure 
measurements are available. 

Consider first the flow for large values of a. The photograph and 
pressure curve for a=2-0in. both confirm the qualitative description 
given in the introduction. In fact, the two patterns at the upstream 
and downstream ends can be treated as independent. The gradual 
decrease in pressure in the intermediate region with attached flow was to 
be expected, as the pressure on the body for very large values of a, or with 
no second cylinder at all, would settle down to a value very nearly equal 
to the free-stream static pressure, i.e., P=0. 

The value of the pressure coefficient at the base of the first cylinder for 
large values of a (see for example the curve for a=2:0 in.) was found to be 
—0-16. This value of the base pressure coefficient was found by Chapman 
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(1950) to occur on a cone-cylinder at M,=2-0 for a Reynolds number, 
based on body length upstream of the base, of about 10. The Reynolds 
number based on body length in the experiments described here was 
3-4 10°, which in Chapman’s experiments gave a pressure coefficient of 
about —0:-10. However, from fig. 19 in Chapman’s paper it can be seen 
thatthe base pressure coefficient falls very rapidly with increasing Reynolds 
number in this range, and it is not surprising that the Reynolds number 
based on the length of a bluff body should correspond to a considerably 
higher value for a cone-cylinder. The boundary layer on the flat-nosed 
cylinder starts growing on the nose, and its thickness is further increased 
by the laminar separation and reattachment. A comparison based on 
Reynolds numbers with the displacement thickness just upstream of the 
base as typical length (if this thickness could have been determined) 
would therefore probably have given much closer agreement between the 
present results and those obtained by Chapman. 

A comparison of the cone angle of the separated region upstream of the 
second cylinder with the photographs by Mair (1952) of separation of 
laminar and turbulent boundary layers on probes in front of flat-nosed 
bodies, indicates that the boundary layer was turbulent in the present 
experiments. In fact, as discussed below, transition was found to occur at a 
distance of about 0-2 in.—0-3 in. downstream of the upstream end of the gap. 

As a was decreased from a=2:0 in. to the critical value of 1-25 in. the two 
separated regions approached each other, and the shock waves at the end 
of the first and the beginning of the second separated regions would 
coincide at a=1-25in. It will be seen from the pressure curve for 
a=1:3 in. that the second cylinder began to affect the pressure at the base 
of the first cylinder before the critical value of a was reached. 

To discuss the flow pattern in the cases with only one dead-air region it 
is convenient to start with small values of a and follow what happened when 
@ was increased towards the critical value. 

For a—0-l in. and 0-2 in. the pressure was found to be only slightly 
different from the free-stream static pressure. On the photographs (see 
fig. 4 (Plate 1)) there were weak compression waves at the beginning and 
the end of the mixing region, showing that this region was slightly 
convex outwards. This could not be concluded from the pressure curves. 

For a between 0°3 in. and 0:8 in. the pressure at the upstream end of the 
gap was still only slightly different from p,, but there was a steep rise in P 
at the downstream end. Furthermore, over most of the gap P decreased 
almost linearly (see for example the curves for a=0-6 in. and 0-8 in.). 
The steep rise in the pressure near the end of the gap and the high value 
found at the hole farthest downstream can be explained by the impact of 
the mixing region where it meets the sharp corner at the end of the gap. 
Part of it is deflected outwards past the second cylinder giving rise to a 
compression wave which could be seen springing from the corner in som e of 
the photographs not reproduced here. The other part of the mixing region 
is deflected towards the axis as a region of high velocity along the upstream 
face of the second cylinder thus causing the high pressure in the corner. 
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The pressure curves show that, if the pressure is assumed not to vary 
much in planes perpendicular to the axis (except very near the down- 
stream end of the gap as discussed above), the mixing region must be 
convex outwards. This is confirmed by the fact that for a between 
0:3 in. and 0-8in. the photographs showed compression waves both at 
the beginning and at the end of the mixing region. 

When a was increased above 0:8 in. the shape of the pressure curves 
changed, and for a=1-2 in. P increased almost linearly over most of the 
length of the gap. This means that the mixing region was now concave 
outwards. This conclusion about its shape is in agreement with the 
photographs. On fig. 6 (Plate 1) where a=1-2 in. the expansion at the 
first corner (which is also evident from the pressure curve) and the com- 
pression wavelets along the mixing region are visible. The concentrated 
compression at the downstream end of the gap can also be seen. It was 
followed by a (hardly visible) expansion bending the flow parallel to the axis. 

The fact that the mixing region was curved outwards and P decreased 
with increasing 2 for small values of a, whereas the mixing region was 
curved inwards and P increased with x for larger values of a, may perhaps 
be explained by the transition from laminar to turbulent flow in the mixing 
region, which could be seen on some of the photographs to take place about 
0-2 in.—0-3 in. downstream of the first cylinder. 

An increase in pressure in the direction of the flow along the surface of 
the gap, i.e. from right to left, must be an inertial effect and correspond to a 
decrease in velocity. On the other hand, a decrease in pressure in the 
direction of the secondary flow may be an inertial effect and correspond 
to an increase in velocity, or it may be due to resistance caused by friction 
and turbulence. For small values of a the flow was probably not very 
turbulent and the resistance was only of moderate importance, whereas 
for larger values of a the pressure distribution was governed by the 
tendency of the turbulent losses to be balanced by a pressure gradient. 

It will be seen from the above discussion that there was a gradual 
change in the flow pattern as a was increased towards its critical value. 

Although the flow in the dead-air region is not known in detail, an 
attempt has been made in fig. 17 to sketch a possible flow pattern for the 
case of a=1-2in. It should be emphasized that this sketch is largely 
based on conjectures and is only intended to show qualitatively how some 
of the observed phenomena may be explained. It shows the mixing 
region curved inwards and drawing air from the dead-air region. At the 
end it is split up into two parts, one passing over the body, the other 
being added to the dead-air region. 

The bump in the pressure curve for a=1-2 in. with a local maximum in 
P at x=0-1 in. is explained by a small secondary vortex in the corner. 
This would give a stagnation point with a relatively high pressure at 
Cas Ue an 

Although no actual values of the velocities in the dead-air region can 
be deduced from the experiments, the great changes in pressure along the 
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gap show that they may be of considerable magnitude. The difficulty 
in giving absolute values of the velocities arises from the fact that the 
effects of viscous and turbulent shear stresses are unknown. Some of the 
photographs indicate that the flow was highly turbulent, and the Reynolds 
stresses may well play a very important part in the production of the flow 
pattern. It would therefore be very desirable to make measurements in 
the dead-air region with hot wires, but this was not possible with the 
experimental set-up described here. . 

In order to determine the actual values of the drag caused by the gaps 
it would be necessary to measure the pressure distributions on the up- 
stream and downstream faces of the gaps. It may, however, be concluded 
from the experiments, that on the body with d=} in. the drag was small 
for a smaller than 0-3in. It then increased but stayed nearly constant 
for a between 0-4 in. and 1-0in. It then again increased as a was increased 
to the critical value of 1-25in. Above the critical value the drag was 
nearly constant independent of the length of the gap. 


Fig. 17 


Sketch of flow for a=1:2 in. 


Tn conclusion it is interesting to compare the results described here with 
some experiments by Kavanau (1954) on the flow at the base of sting- 
supported bodies of revolution at a Mach number of 2:84 and low Reynolds 
numbers. Kavanau used a supporting sting which was continued after a 
certain length in a much thicker cone-cylinder. — The critical sting length 
was defined as the minimum sting length possible for obtaining a base 
pressure within 0-5% of that obtained with an infinitely long sting. 
This length is therefore slightly greater than the critical value of a as 

i resent paper. 
Se ain ine the critical length increased rapidly with decreasing 
Reynolds number based on body length. He also plotted the base 
pressure against sting length. A similar curve would be obtained by 
plotting P forz=0 asa function of a as determined in these experiments. 
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SUMMARY 


For an infinitely-long circular cylinder oscillating rectilinearly along 
a diameter, the damping force is determined from experimental measure- 
ments of the damping force on a circular cylinder oscillating about one 
end. At very small amplitudes of vibration, up to about one-tenth of 
the diameter, good agreement is obtained with Stokes’ theoretical result 
for the viscous damping force at infinitesimal amplitudes of vibration. 
But at larger amplitudes, the damping force is found to increase linearly 
with amplitude ; an empirical theory is given which is in qualitative 
agreement with this result. 


$1. LyrropucTION 


THE viscous damping force on an infinitely-long cylinder vibrating 
rectilinearly normal to its length with infinitesimal amplitude has been 
calculated by Stokes (1851), who showed that the force is proportional 
to the product of angular frequency w, velocity U and mass of fluid WM’ 
displaced by the cylinder. A boundary-layer treatment of the problem 
by Stuart (1952) gives a value of the coefficient of proportionality k’ 
which is half that of Stokes’ formula. Martin (1925) has obtained values 
of the damping on vibrating strings of circular cross section, and his 
results agree well with Stokes’ result. His experiments were carried 
out with small amplitudes of vibration, and he states that no influence 
of amplitude on the damping could be detected. On the other hand, 
measurements on bridge models by Scruton (1952) show a very marked 
increase of damping with amplitude. Accordingly the experiments 
described in this report were devised to test the validity of the theoretical 
formulae and to determine the influence of amplitude. 

For simplicity, the damping was measured on a long circular cylinder 
oscillating as a pendulum about one end. For amplitudes at the lower 
end of the pendulum of up to about one-tenth of the diameter, no effect 
of amplitude on the aerodynamic damping could be detected, and the 
results are in good agreement with Stokes’ formula. Martin’s (1925) 
amplitude of vibration was of the same order of magnitude. In the 
present tests, it was found that at very much larger amplitudes than this, 


the pressure damping far exceeds the viscous damping. 
ee Ea eee 
* Communicated by W. P. Jones. 
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Additionally, the result for the viscous damping at small amplitudes 
confirms the view that the boundary-layer method used by Stuart (1952) 
is not valid for the calculation of viscous damping, because the pressure 
distribution is modified by viscosity by an amount of the same order of 
magnitude as the damping force. 


§ 2. DESCRIPTION OF APPARATUS AND MretrHop oF MEASUREMENT 
2.1. Measurements at Low Amplitudes 


The apparatus consisted of a circular cylinder (with a length to diameter 
ratio of 24) suspended at one end by a cross-spring hinge from a metal 
plate, the whole being mounted in a cylindrical pressure tank to enable 
experiments to be made at several pressures. In addition, the pendulum 
could be suspended in free air, in order to examine the aerodynamic 
effects of the pressure tank walls. 

The movement of the pendulum was indicated electrically on a double- 

beam cathode-ray oscillograph, and the trace together with a suitable 
time base was recorded on moving film. Five records were taken, one in 
free air at atmospheric pressure, and four in the tank at various pressures. 
From each record the corresponding logarithmic decrement was found, 
and this was constant within the limits of experimental accuracy for 
amplitudes less than 0-004 radians. The frequency f also was measured 
from the film records, while the moment of inertia J was determined from 
stiffness and frequency measurements. The aerodynamic decrement 
5, was obtained by subtraction from 6 of the extrapolated (apparatus) 
damping 5, at zero pressure. Results are given in table 1. 


Table 1 


NNN. 


Pressure : 
Location genet Frequency a Log. dee. Aerodynamic 
f ; : i : log. dec. 
: inertia inches of ) 
cylinder T Ib ft? C.p.s. mercury s 
Free air | 10-60 | 0-690 | 30-30 | 0-00736 | 0-00223 
10-60 0-700 58:70 | 0-00842 0-00329 
k 10-60 0-702 30-59 0-00722 0-00209 
aca 10-60 0-703 14-66 0-00646 0-00133 
10-60 0-703 0-54 0-00520 — 


‘7 en ne nnn 
2.2. Measurements at Large Amplitudes 


A spot of light was projected onto a screen on the floor via a mirror 
on the pendulum near to the pivot. An average logarithmic decrement 
5, for the amplitude range was then obtained from the time ¢ paeen fop 
the oscillation to decay between two given amplitudes 6 and @ » and 
from the frequency f; the aerodynamic decrement by was given by 
subtraction of the known apparatus decrement 8), which was assumed 
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to be independent of amplitude. The experiment was repeated for 
several different average amplitudes. 


0-010 
0-008 
dv 


0-006 


0-004 


Aerodynamic logarithmic decrement against angular amplitude for large 
amplitudes. 


Results are given in table 2, while fig. 1 shows the quantity 5, plotted 
against the average amplitude 6. The curve is drawn to pass through 
the logarithmic decrement 5, at small amplitudes. An empirical theory 
for large amplitudes is given in an Appendix. 


Table 2 
Amplitude range Mean amplitude : : f 
O=(6)+0,)/2 if 
rads. rads. rads. C.p.s. 
0-079 0-060 0-070 0-0171 0-0119 0-692 
0-059 | 0-040 0-050 0-0130 0-0079 0-692 


0-040 0-020 0-030 0-0113 0-0061 0-692 
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§3. COMPARISON WITH THEORY FOR SMALL AMPLITUDES 


It is first necessary to relate the logarithmic decrement of the pendulum 


with Stokes’ damping function for an infinitely long cylinder vibrating 
normal to its length. ; 


Fig. 2 


Notation. 


The notation for the circular cylinder is shown in fig. 2. Consider an 
element of length dx, and assume that the aerodynamic force on the element 
is given by k’wUM’ (see § 1) for an infinite cylinder vibrating normal to 
its length. Then the damping force is mparxbw(k'+hk") dx, where k’ is 
the aerodynamic damping coefficient, k" is the apparatus damping 
coefficient, and @ is the angular velocity. The total moment about the 
pivot of the damping force is found by integration to be 


(—4)prra®lcu(h’ +h") 8. 
Thus the equation for the pendulum’s motion is 
T6--tpra%u(k'tk")O+osn0=0, . .. ~~ (I) 
where I is the moment of inertia, including the aerodynamic inertia, 
o is the stiffness, and @ is the angular acceleration. Putting [=4p,7a7/, 
where p, is an equivalent density, and o=Iw?, we have 


b+ = afl E+ ot =0, Swe he) «2 (2) 
1 


for small amplitudes. Thus the total logarithmic decrement 6 is given by 
S=(p/p,)7(k’ +k"), while the aerodynamic logarithmic decrement is 


6 y= ee e . e ” ° ° ° . (3) 
P1 


From experiment § y=0-00223 and. py=22°8 lb tt, while p=0-0765 Ib fe", 


Thus k’=0-211. 
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When the parameter a1/(w/v) is sufficiently large, Stokes (1851, eqn. 
(7.27)) gives the damping function k’ as 


2 
fe a a 


In the present case ay/(w/v) is 13-8, which certainly lies within Stokes’ 
asymptotic range. Substitution of various quantities into (4) gives 
k'=0-215, which is very close to the experimental value of 0-211. But if 
the second term in (4) is neglected and the assumption k’=(C/a)1/(v/w) 
is made, where C is a constant, C may be calculated from the experimental 
k’ to be 2-94, which is in good agreement with the value 2,/2 of Stokes 
(1851) and the value 2-87 obtained experimentally by Martin (1925). 

Since the viscosity, w~=vp does not vary much with density, it is seen 
from (3) and (4) that 6, is theoretically proportional to p/?, where p is 
the air density, or to p’/?, where p is the pressure. The values of the 
logarithmic decrement in the tank (table 1) for different pressures are 
only in qualitative agreement with this. The ratio of 5, at a pressure 
of 58-70 in. of mercury to that at 30-59 is 1-57, and the ratio of 5, at a 
pressure of 30-59 in. of mercury to that at 14-66 is also 1-57, while each 
should approximately be 1/2. This discrepancy may be due to the 
influence of the tank wall. 


$4. ConcLUDING REMARKS 


The range of amplitude for which the damping is effectively viscous 
and independent of amplitude is found to be about 0-004 radians, which 
at the lower end of the pendulum corresponds to one-tenth of the diameter. 
Above this amplitude, pressure damping increases in importance, as 
shown by the rough experiment done to determine the damping at large 
amplitudes (fig. 1). An empirical, quasi-steady theory for large ampli- 
tudes is presented in an Appendix. In spite of the extravagant nature 
of the assumptions made, the theory predicts a linear increase of damping 
with amplitude of the same order as that found by experiment. 

Before reaching the above conclusions about the results for the small 
amplitude tests, the following factors were considered : 

(a) The effect of pressure variations on the results obtained inside the 
tank is not quite in accordance with the square root law. 

(b) The wall of the tank may have an influence on the damping even 
at low amplitudes ; this view is supported by the comparison of the 
logarithmic decrement at atmospheric pressure inside the tank (0-00209) 
with that in free air (0-00223). This may account for the divergence 
noted under (a). But, of course, the tank wall has no effect on the 
extrapolated value of the apparatus damping. 

(c) The transfer of results from the case of a cylinder vibrating about 
one end to the case of a cylinder vibrating normal to its length neglects 
the longitudinal acceleration in the former case. But if 6) is the maximum 
angular amplitude of the pendulum the maximum tangential acceleration 
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of the fluid is of order /w6), while the maximum longitudinal acceleration 
is of order lw6?; the ratio of longitudinal to tangential acceleration 
is thus 0). This is required to be small compared with unity, and in fact 
was about 0-004 radians in the small amplitude tests. 

(d) In view of the large length to diameter ratio of the pendulum, end 
effects seem unimportant. 

The authors are indebted to Mr. C. Scruton for his continued interest 
and helpful suggestions during the course of this work. 

The work described above was carried out in the Aerodynamics Division 
of the National Physical Laboratory, and this paper is published on the 
recommendation of the Aeronautical Research Council and by permission 
of the Director of the Laboratory. 
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APPENDIX 


Empirical Treatment of Damping at Large Amplitudes 

The logarithmic decrement of an oscillation is normally defined as the 
ratio of amplitudes of successive cycles. It can also be expressed as 
an energy decrement. Suppose an oscillation of energy H at a given 
amplitude suffers a loss of energy 4H due to damping during the next 
cycle; then the logarithmic decrement 6 is 4H/2#, provided 6 is small 
compared with unity. 

For the case of an oscillation of large amplitude, it is now assumed 
that the damping force or drag on the cylinder is given by the quasi- 
steady expression CU V? per unit length, where C is a constant and V is 
the local velocity at a point of the pendulum. In the notation of §3, 
the work done by the drag per cycle is 


1 0/200 : 
An=4| | Gace Oy 
2=0/ t=0 


Since 8 is small, we can neglect the amplitude decrement during one cycle 
and write 0=6, sin wt, where Oy 1s the angular amplitude at the particular 
stage of oscillation. It then follows from (A 1) that 4H=(2)Cl*w6,°, 
and since H=(})[w6,?, it follows that 

2cr 


3y= 3-7 ee 4D) 


Since the damping force is assumed to have the quasi-steady form 
CYV2, C is given by steady motion results (Goldstein 1938). In the large 
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amplitude tests summarized in fig. 1, 0) varies between about 0-07 and 
0-03 radians. Thus the Reynolds number based on diameter has a 
maximum value of between 1270 and 540 at the lower end of the pendulum. 
Reference to Goldstein (1938) then shows that the drag coefficient Cp 
may be taken as unity, so that C=ap, where a is the cylinder radius and 
p the air density. 

On substituting in (A 2) the values a=1/12 ft, J=4 ft, [=10-60 lb ft? 
and p=0-0765 lb ft-*, we obtain for the aerodynamic logarithmic decre- 
ment 

S,==0'1020,5°8. tia as 1 ee 
compared with the value 
8,==0:153 05 >.) e a ie) Lauren ceed 
indicated by fig. 1. 

This empirical theory thus shows (eqn. (A 2)) that the logarithmic 
decrement, 6,,, is proportional to the angular amplitude, 65, as indicated 
by experiment to be the case approximately. In addition, the numerical 
- value obtained for 6, is in fair agreement with experiment. 

The assumptions made are :— 

(i) That the logarithmic decrement is small compared with unity, an 
approximation which is justified by the magnitude (0-01) of the experi- 
mental values of 5,. For far higher amplitudes, the expression 6 ,= 42/2H 
requires modification. 

(ii) That the damping force is given by a quasi-steady expression, 
which entails neglect of acceleration. 

(iii) That a single constant, C, which in the present case is based on a 
maximum velocity, is sufficient to specify the drag in the form CV? 
at all velocities. 


Lona od 
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ABSTRACT 


A comparison is made between the heat conductivity data below 1°K in 
which the specimens were exposed to the magnetizing field and those 
observations in which the specimens were shielded. The thermal 
resistance of a tin single crystal in which superconductivity was destroyed 
by a transverse magnetic field has been examined in detail. The results 
are discussed and an explanation for the heat conductivity minimum in 
the intermediate state is suggested. 


$1. IytTRODUCTION 


Tue heat conductivity of pure superconductive metals has in recent years 
been the subject of a great deal of investigation and a fairly clear pattern 
of its behaviour has been evolved. K,,, the heat conductivity in the 
normal state, measured below the critical temperature 7’, by destroying 
superconductivity with a magnetic field, is found to pass through 7’, 
without a break. K,, the heat conductivity in the superconductive state 
on the other hand is found to depart from K,, at T,.. W,, and W, which will 
be referred to are the corresponding heat resistivities. 

When a transverse magnetic field is applied to a superconducting 
cylinder, penetration commences when the applied field H is equal to 
1/2H,, and is complete when H=H,; H, is the critical field. Between 
H=1/2H, and H, the specimen is said to be in the intermediate state. 
The exact structure of this intermediate state has been the subject of 
many theories and investigations which in a rather over-simplified form 
suitable for use here reduces to a pattern of alternate layers of normal and 
superconductive material. Then, if the fraction of the length which is in 
the normal state is proportional to the strength of the penetrating field, it is 
easily seen that the thermal resistance W should change linearly from 
W, at H=1/2H, to W,, at H=H,. This behaviour was indeed found in 
the first experiments (Mendelssohn and Pontius 1937, de Haas and 
Rademakers 1940) but a quite unexpected effect was found by Mendelssohn 
and Olsen (1950 a, b) in specimens of columbium and lead with a trace of 
ee esi 

* Communicated by Dr. K. Mendelssohn, F.R.S. 
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bismuth. In these, when entering the intermediate state, W first increased 
to a maximum before falling to W,. This has more recently been 
observed for pure metals (Webber and Spohr 1951, Detwiler and Fair- 
bank 1952). 

The experimental investigations which have been made by us have been 
chiefly concerned with the temperature variation of K,, for some pure super- 
conductors below 1°K. This has been fully reported elsewhere (Men- 
delssohn and Renton 1954) and it was observed that distinct differences in 
the experimental results are found, depending on whether or not the speci- 
men is shielded from the applied magnetic field used for the cooling process. 
The results of the preliminary work not employing magnetic shielding were 
reported earlier (Olsen and Renton 1952, Mendelssohn and Renton 1953). 
The explanation of this discrepancy lies in the fact that there may be some 
‘frozen in’ field in the specimen after demagnetization and the resulting 
measurements will then refer to the heat conductivity in the intermediate 
state. 

The thermal resistivity maximum in the intermediate state was found 
in a lead single crystal (Olsen and Renton 1952) at three temperatures. 
The specimen here was not shielded and the results are inconclusive, hence 
a further experiment has been performed using a shielded tin single 
crystal. 

In this paper a comparison is made between these preliminary results 
and the results obtained when the specimen is magnetically shielded. We 
also give the results of a field transition cycle and in the discussion we 
propose a possible explanation of the occurrence of this maximum of the 
thermal resistance in the intermediate state. 


§2. REsuLts 


Figure 1 shows a transition cycle between the superconducting and 
normal state for the tin single crystal specimen the heat conductivity of 
which had been measured earlier. The specimen which is shielded from 
the demagnetizing field starts with its heat conductivity in agreement with 
the results given in Mendelssohn and Renton (1954), exhibits a maximum 
in W in the intermediate state, and on returning to zero field shows 
evidence of a ‘ frozen-in ° field which reduces its conductivity by a factor 
two in agreement with the preliminary result of Mendelssohn and Renton 
(1953). On reversing the field the repeated cycle no longer shows the 
maximum in the intermediate state. 

The general pattern of the comparative results for K, for the various 
specimens with and without shielding is as follows. 

Lead Single Crystal. The heat conductivity is the same in both cases. 

Tin Single Crystal. The heat conductivity is higher with shielding than 
without. 

Tin Polycrystal. The heat conductivity is higher with shielding than 
without. 
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The change in the thermal resistance W is shown as superconductivity is 
destroyed by a transverse field. The specimen used was the tin single 
crystal whose conductivity had previously been measured below 1°K. 
This cycle was performed at a temperature of 0-34°K where K, is purely 
lattice conduction. 

Fig. 2 
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Indium Single Crystal. The heat conductivity of this specimen is 
higher than that of a polycrystalline specimen of the same source, both of 
which were magnetically shielded. Both these specimens in turn had a 
higher conductivity than another single crystal made from metal of the 
same source which was not magnetically shielded. 

Thallium Polycrystal. Of all the specimens investigated, here alone the 
heat conductivity for the shielded specimen was lower. It should be 
noted, however, that the specimen used in the preliminary work had been 
destroyed by oxidization. The second specimen which was shielded was 
made from metal of the same source. 

Tantalum Polycrystal. The heat conductivity is higher with shielding 
than without. The results for this specimen are shown on a double 
logarithmic plot in fig. 2. The results for the shielded specimen are 
derived from three separate experiments whilst the two unshielded results, 
showing different amounts of ‘ frozen-in ’ field, are each the result of one 
particular cooling. 

Columbium Polycrystal. The heat conductivity is higher with shielding 
than without. 


§ 3. Discussion 


The results for the specimens which have not been magnetically shielded 
fall readily into place if it is considered that the specimens contain a 
frozen-in field. The measurements then refer to a particular stage of the 
intermediate state such as values lying between H/H,—0-5 and 1 in fig. 1. 
Since in all cases except lead and thallium the conductivity under these 
conditions is reduced, it appears that these specimens would all exhibit 
a maximum in the thermal resistance in the transition from the super- 
conducting to the normal state. In the case of lead alone it appears that 
there is no frozen-in field, whilst, as has already been discussed, the two 
specimens of thallium are not really comparable. 

No explanations of the maximum in the intermediate state can be made 
by any additive process. The explanations which have been put forward 
have been comparatively discussed by Detwiler and Fairbank (1952) who 
come to the conclusion that the most probable explanation lies in the 
assumption of some additional scattering mechanism of either lattice waves 
or electrons at the superconducting—normal boundaries of the laminae. 
Two further accounts have since appeared due to Cornish and Olsen (1953) 
and to Hulm (1953) both of which are open to slight objections. Cornish 
and Olsen propose a model of a lattice and electronic system side by side 
having effectively a different temperature distribution and with an 
interaction proportional to the temperature difference between them. 
Such a detailed separation into two distinct systems may be over- 
simplified (for instance see Mendelssohn (1953), discussion) and can really 
only be used where the conduction due to one system is of a different order 
of magnitude to the other. 
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at the additional thermal resistance is due to either scattering of electrons 
or lattice waves at the superconducting—-normal boundaries, then the 
amount by which the heat conductivity is decreased will give a minimum 
; estimate for whichever conduction mechanism is considered. Now, 
electronic conduction predominates at the highest temperatures and 
lattice conduction at the lowest temperatures and since it is found that the 
effect increases at lower temperatures, and in most cases does not occur at 
temperatures near 7’, at all, this would seem to point to the effect being 
due to the lattice component. Moreover, the effect can be observed at the 
lowest temperatures where conduction is known to be almost exclusively 
lattice conduction (Mendelssohn and Renton 1954, loc. cit.). In this 
temperature range it would be impossible to account for the increase in 
resistance as due to the scattering of electrons at these boundaries. This 
collected evidence would favour the explanation in terms of a lattice 
resistance at the superconducting-normal boundary and seems therefore 
to rule out such attempts as have been made (Webber and Spohr 1952, 
Hulm 1953) to explain the maximum as a decrease in electronic conduction. 

The way in which such an increased lattice resistance will occur can 
be seen from the Casimir (1938) theory of the boundary scattering of 
lattice waves. Here the wavelength of the lattice waves is limited by the 
physical boundaries of the specimen. If, as will be discussed below, the 
superconducting—normal boundaries act in the same manner as physical 
boundaries, the effect can easily be understood. In the intermediate state 
an increase in the overall length of the normal regions tends to decrease the 
thermal resistance of the specimen whilst an increase in the number of 
superconducting-normal boundaries tends to increase the resistance by 
reducing the mean free path of the lattice waves. The progression into 
the intermediate state involves the relative rates of change in these oppos- 
ing mechanisms. At first the number of boundaries increases to increase 
the thermal resistance more rapidly than the change in length of the normal 
regions decreases the resistance. The resistance rises to a Maximum and 
then falls when the overall length of the normal regions changes more 
rapidly. After this maximum the decrease in the resistance probably 
changes more rapidly still if the normal regions expand into one another, 
eliminating some of the superconducting laminae rather than forming new 
normal laminae within existing superconducting laminae. 

To allow for this explanation of the maximum, it is necessary to sub- 
stantiate our postulate that a superconducting—normal boundary acts in 
the same way as a physical boundary. At first sight this would seem 
unlikely since the difference between the two states is purely electronic. 
As recent experiments (Mendelssohn and Renton 1954) have however 
shown, there exists an important indirect effect due to the much reduced 
scattering of lattice waves by electrons in the superconductive state. 
On the other hand, in the normal state scatter by electrons was found to be 
the predominant process in limiting lattice conduction. Thus, in the 


normal regions of the intermediate state where the density of normal 
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electrons is high, lattice waves will be subject to a great deal of scattering. 
In the superconducting regions where the density of normal electrons is 
low and gets lower with decrease of temperature, the lattice waves will 
hardly be scattered. For the lattice waves the normal regions will there- 
fore act as a dense medium while the superconducting regions will be 
transparent. The degree of this transparency increases as the temperature 
falls and in this way the boundary between the two media can limit the 
mean free path of the lattice waves. It may be that the magnitude of the 
effect increases at lower temperatures not only because the proportion of 
lattice conduction increases but also possibly because the transparency of 
the superconducting region relative to the normal region increases. This 
would then make the boundaries a more effective limit to the lattice waves. 

The relevance of the comparison of the results with and without a 
frozen-in field is of interest in considering the superconducting heat switch. 
It is seen that in most cases a frozen-in field reduces the heat conductivity 
of the specimen and that most specimens show evidence of such a frozen-in 
field. In the use of a heat switch the superconductor is first rendered 
normal with a magnetic field and then superconducting by removing 
this field and so the conditions for producing such a frozen-in field occur. 
Thus, conveniently, the efficiency of this switch is greater than expected 
since it is required that K, should be as low as possible. It has been found 
in practice (Darby, Hatton, Rollin, Seymour and Silsbee 1951) using a 
lead wire that K, was indeed a factor 60 smaller than the value extra- 
polated from results at 2°k. This could be due to both a frozen-in field 
and also due to the fact that whilst the value at 2°K is predominantly 
electronic conduction which is known to fall off rapidly, below 1°K the 
conductivity will be limited by lattice conduction. Thus below 1°K the 
conductivity will be limited by the size of the crystallites of the specimen 
and any frozen-in superconducting—normal boundaries which further 
reduce the size of the internal structure. 
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VII. Experiments on Supersonic Flow Over Flat-Nosed Circular Cylinders 
at Yaw—II: Pressure Measurements on a Cylinder at 10° Yaw 
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SUMMARY 


In a previous paper (Hall 1954) the supersonic flow over flat-nosed 
circular cylinders at yaw was investigated by schlieren photography. 
This investigation was continued for the case of a cylinder at 10° yaw by 
measurements of the surface pressure distribution, which are described 
in the present paper. 

Close to the nose the cross force, obtained by integrating the pressure 
distribution, was found to be large. This was due to the asymmetry of 
the boundary layer separation. Further downstream the variation of 
the cross force with distance from the nose was in good qualitative agree- 
ment with the variation with time of the drag force on a two-dimensional 
circular cylinder started impulsively from rest. 


§ 1. IyrrRopuUcTION 


In an earlier paper (Hall 1954) the flow pattern about a flat-nosed circular 
cylinder at yaw to a supersonic stream was deduced from schlieren 
photographs. An analogy was drawn between the development of the 
cross flow with distance from the nose and the development of the two- 
dimensional flow pattern around a circular cylinder started impulsively 
from rest. 

Additional information has since been obtained from pressure measure- 
ments on the surface of a cylinder at 10° yaw. The results are described 
below. 

§ 2. EXPERIMENTS 

The experiments were made in the same wind tunnel and under the 
same conditions as in the first paper. 

The model was a straight steel tube 5/16 in. outside diameter with a 
thin disc soldered on the upstream end. Pressure holes were drilled at 
points on a generator of the cylindrical surface and on a diameter of the 
nose. The centres were all in one plane. The holes were of 0-015 in. 
diameter except for one hole near the edge of the nose which was of 0-005 
in. diameter. They were filled with solder and each hole was redrilled. 


as required. 


A a a aS 
* Communicated by P. R. Owen. 
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The model was supported in the wind tunnel with its axis inclined at 
10° to the free-stream direction. By rotating the model about its axis 
the pressure holes were set at different values of , the angle between 
the plane containing the centres of the pressure holes and the plane of 
symmetry of the flow.* The pressure, p, was measured at each hole for 
several values of ys. 

The pressure distributions on the nose and on the cylindrical surface 
are presented separately using two different non-dimensional pressure 
coefficients. This facilitates comparison with work by other authors. 
The pressure distribution on the nose is given in terms of p/H, where H, 
is the stagnation pressure behind a normal shock wave in the free-stream 
of Mach number 1-965, calculated from the measured stagnation pressure 
upstream of the throat, H. 


Fig. 1 


Pressure contours on the nose. Values given are of p/H,. 


For the distribution on the cylindrical surface the non-dimensional 
coefficient C,=(p—p,)/4pU is used. p,, p, and U are the static pressure, 
the density, and the velocity in the free stream. At a free-stream Mach 
number of 1-965 the pressure coefficients are related by 

P 


C,=2-020 = —0-370. 
(p= 2-020 FF — 0-370 


§ 3. Tor PressuRE DisTRIBUTION ON THE NosE 


The pressure distribution on the nose of the cylinder at 10° yaw is 
given in fig. 1. The position of the stagnation point can be determined 
approximately from the pressure distribution along the diameter in the 
—_—_—- 

* b=0° is on the pressure side of the body. This is not consistent with the 
previous paper but the change was made to conform with the usual practice. 
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plane of symmetry (fig. 2). It occurs at about 0-17d from the axis, 
where d is the diameter of the cylinder. At the stagnation point 

1g <—_ . 

He == 0-984.* 

Since H, is the minimum stagnation pressure just behind a shock wave 
in a stream of Mach number 1-965, it is evident that there must have been 
a loss of total pressure somewhere upstream of the body, in addition to. 
the known loss through the shock wave. A similar loss was observed 
by Murphy (1953) during the calibration of supersonic wind tunnels with 
a pitot tube. It is hoped to make further investigations to try to explain: 
the mechanism of this loss. 


Fig. 2 


0.4 0.2 0 02 0.4 
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p=180<— 7 — >= 0 


The pressure distribution on the diameter of the nose in the plane 
of symmetry of the flow. 


The coefficient of normal force on the nose may be defined as 

_- Force—jnd? p, 

Or 5 U2 ind 

On p y) 
= y dy db— a 

J0 He Gn yM* 


(72) 
to 
—_ 
——_— 
i) be 
= 
Ww 
See 
bo 


where y is the distance in diameters from the centre of the nose. The 
integral was evaluated with the assumption that p/H, Cee! to the 
sonic value of 0-528 at the edge of the nose (Holder and Chinneck 1954). 
This value may be in error but the precise value is not important for the 
evaluation of the integral. It was found that 


Cp Z 1°38. 


* The last figure should not be taken as significant but the value of p/H, may 
be assumed with certainty to lie between 0-98 and 0-99. 
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The value of Cy, when the body is at zero yaw was determined by extra- 
polating the results of Holder and Chinneck, in the Mach number range 
1-4 to 1:8, to a Mach number of 1-965. In this case 


Cp & 1-43. 


§ 4. Tok PRessURE DISTRIBUTION ON THE CYLINDRICAL SURFACE 

A point on the cylindrical surface of the body is specified by the 
coordinates x, %, where x is the distance from the centre of the nose 
measured along the axis of the body. The distribution of pressure with 
% was measured at several values of x. The results are given in figs. 3 
to 7. 

The pressure distribution close to the nose is governed by the boundary 
layer separation at the shoulder and the reattachment further downstream. 
Inside the separation ‘ bubble ’ there is a ° dead-air ’ region in which the 
pressure is approximately constant. This is shown in figs. 3 and 4, 
which show pressure curves at stations close to the nose. At x’ =a/d=0-06 
the pressure is almost constant for all % indicating that the boundary 
layer is separated all round the body. Further downstream the constant 
pressure region does not extend over all %, showing that reattachment 
has taken place for some values of #. 

On the pressure side of the body the pressure increases after the 
reattachment of the boundary layer. The pressure curves for 2’ =constant 
and for % between 0° and 90° tend towards a curve which is in good 
agreement with the pressure curve calculated with the assumption of a 
potential cross flow. The theoretical curve is shown as a dashed line 
in fig. 4 and represents the following pressure distribution : 

C,,=(1—4 sin? #) sin? 10°. 

Between a’=1-18 (fig. 4) and x’=2-46 (fig. 5) there is a large increase 
in pressure near 4—180°. This was observed in the schlieren photo- 
graphs as a shock wave (Hall 1954). On this side of the body the region 
adjacent to the body is essentially a * dead-air’ region for x less than 
1-18, whereas, for x’ greater than 2-46 it is essentially a vortex wake. 
In the latter region, the wake consists of two symmetrical vortices which 
induce a flow towards the body. This results in a high pressure near 
#=180°. For x’ greater than 2-78 the pressure near %=180° falls 
slowly with increasing w’. This is attributed to a decrease in the maximum 
flow velocity between the vortices as they become further apart. Eventu- 
ally, the pressure peak disappears completely («’=9-50, fig. 7) where 
the wake no longer consists of two steady, symmetrical vortices. 

The wake is not formed in exactly the same manner as would be 
expected from the two-dimensional analogy. In the latter case the 
initial separation occurs at the rear stagnation point. The separation 
point moves forward and the vortices are formed. In the case of the 
yawed cylinder the separated boundary layer from the nose does not 


reattach near = 180° and merges into the beginning of the wake. This 
is illustrated in fig. 8. 
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Figs. 3 and 4 


Pressure distributions at stations on the curved surtaze. 
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Figs. 5-7 


oe 45° 902 135° 180°" 


Pressure distributions at stations on the curved surface. 


This coefficient is analogous to the corresponding dra 
two-dimensional flow past a circular cylinder. 
The variation of Cgp with x’ is shown in fig. 9. 
nose is caused by the asymmetry of the separation bubble. 
boundary layer is attached (on the pressure side) 
whereas there is a much lower pressure where the la 
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The pressure distribution in the analogous two-dimensional flow about 
a cylinder started impulsively from rest has been obtained by Schwabe 
(1935): the results are given by Goldstein (1938). His pressure curves 


closely resemble those in figs. 5 and 6. In addition to showing the pressure 
rise near 4=180° they also exhibit similar undulations between 4=90° 
and J=144°. 


emt 


A force coefficient, Cop due to the cross flow drag, is defined by 


| (p—p,) cos pdx 
Cop= — 


$oU? sin? 10° 


Fig. 8 
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Flow phenomena adjacent to the curved surface. 


g coefficient for the 


The peak near the 
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there is a high pressure 
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downstream the force coefficient increases to approximately 
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falls again. A similar peak value (2:07) was found by Schwabe in the 
two-dimensional case. 

Unfortunately, it was not possible to make measurements far enough 
downstream for the cross flow to have become steady. It is likely that 
Cop would fall to approximately 0-35 since the boundary layer would 
probably be turbulent at this station. The last point on the curve is 
consistent with this decrease. 


Fig. 9 


x /d ——> 


The variation of cross force with distance from the nose. 
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VIII. Neutron Deficient Isotopes of Pb and TI—I : Isomers of Pb and Ti 
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From the point of view of isomerism all the isotopes of Pb and TI are 
interesting. The odd isotopes of Pb should have isomeric transitions of 
the type i43/5 > f5;. (N<126). 2°6Pb, 2°¢Pb and 2°2Pb have all been shown 
to have high spin isomeric states (Alburger and Pryce 1954, Maeder ei al. 
1954, Hollander et al. 1953) and the same may be true for even—-even Pb 
isotopes of lower masses. Tl has 81 protons and its odd isotopes might 
have isomeric transitions of the same type earlier found in Au (Gillon et al. 
1954) (Z=79). Recently a 0-548 activity has been reported in Tl 
(Henrikson ef al. 1953), which may belong to this class. In the odd—odd ~ 
isotopes of Tl both the odd neutron and the odd proton can be in high spin 
states making a resulting high spin state possible. Recently such a case 
was found in !°°™T'] (Bergstrom et al. 1953). 

For this reason an extensive B-spectrometer investigation was started at 
the Nobel Institute of Physics, Stockholm and at the Gustaf Werner 
Institute for Nuclear Chemistry, Uppsala. Though not yet completed it is. 
worth while to point out some of the most striking features of this 
investigation, reported more in detail in the two following letters. It 
should be emphasized that the recent calculations of the L-conversion 
coefficients by Rose ef al. highly facilitated the interpretations of the 
results. 

Pb. Even—odd isotopes.—There is only one i43/.—>f,,, transition known 
with certainty in lead (Hollander et al. 1953). (07mPb H,—1-063 and 
0-565 Mev, t=0-9s.) No activity due to ?™Pb is reported. A 56s 
activity has been assigned to 22™Pb (#,=0-89 mev) and a 50s activity to 
20ImPh (Hopkins 1952). (#,=0-25, 0-42 and 0-67 Mev.) After the 
discovery (Maeder et al. 1954) of 3-5 h ??™Pb the mass assignment of the 
5-6 s activity seems doubtful and it could therefore be due to oP peor 
this reason lead was irradiated with 60 Mev protons in the Uppsala 
synchro-cyclotron. The Bi activity produced was absorbed in a Dowex 
1 column. The lead activity washed through was directly counted with a. 
scintillation spectrometer. In addition to the well-known Lh 2°4mPh 
activity there was a strong indication of an activity of 10s (£,0-9 Mev). 
Its mother substance decayed with a half-life of 10 h and is probably 
203Rj. 204Bi has a similar half-life (12 h) but it is not probable that this is 
the mother substance of the 108 activity, since 204Pb already has an 
isomeric state (1-1 h). The 10s activity should therefore belong to er | 
It may be identical with the 5-68 activity earlier assigned to ? m™Pb. 
Indications were also found for a 0-5-+0-1 Mev y-ray in the 10s activity. 


* On leave from the Instituut voor Kernphysisch Onderzoek, Amsterdam. 
+ Communicated by Professor Siegbahn and Professor Svedberg. 
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If this y-ray could be confirmed it should be a strong indication that the 
ground state of 2°3Pb has a spin 1/2 or 3/2 instead of f;/2 (Prescott 1954). 

The Uppsala experiments (cf. III) showed that 1°°Pb has a 12 min 
isomeric state from which a 423 kev M4 y-ray is emitted. No second 
y-transition has yet been found. There were also some indications that 
the 234 kev y-ray in 197Pb (71 h) may be due to an isomeric transition. 

Though the information still is incomplete the discussion above shows 
that also in the lead isotopes the energy of the i,3/,—>f;;. transition 
increases (decreasing half-lives) towards the magic number N=126 (see 
also Mihelich and de Shalit 1954). 

Pb. Even-even isotopes.—°*Pb, 2°4Pb and ?°Pb have been suggested to 
have 7-, 7~, and 9~ states respectively (Alburger and Pryce 1954, Maeder 
et al.* 1954). The reinvestigation of 2°™Pb mainly confirmed the results 
of Maeder et al. . It has earlier been suggested that 1°*Pb should have a 
25 min high spin isomeric state (Passel et al. 1954). This assumption 
was based on the fact that a Tl activity of }2h assigned to 1°5Tl was 
milked from a lead activity of +25 m (Neuman and Perlman 1950). 
Since a 18h 198™T] activity later was shown to be due to a high spin 
isomeric state the same should be true for the 25min Pb activity. 
However, in our experiments there was no evidence for any isomeric 
198Pb and 2°°Pb activities. If they exist at all they must be short lived. 
It is therefore concluded that the 2 h Tl and the 25 min Pb earlier assigned 
to 198Pb are identical with the 2-8 h 19’Tl and 42 min 1!9’Pb (ef. III) not 
known when the conclusion above was made. 

Recently (Krohn and Raboy 1954) a third y-ray has been reported in the 
‘decay of 2°mPb with an energy 15--5 kev less than that of the isomeric 
transition and the decay scheme given in fig. 1 has been suggested. Recent 
measurements by us confirm the existence of a y-ray with an energy 
13-2-L0-5kev lower than the isomeric 913 kev y-ray (cf. 11); ifit is in cascade 
with the isomeric ray as suggested by the experiments of Krohn and 
Raboy, the intensity of its K-line indicates it to be an E2 transition. 
‘This suggests the spin assignment given in fig. 1. Thus the first two excited 
states are 2+ and 4+ in agreement with the empirical systematics of 
even—even nuclei. It has earlier been found (Maeder et al. 1954) that the 
most reasonable nuclear shell model assignment for an isomeric level in 
204Pb should be 9~ and its excitation energy comparable with the 9~ level 
in Pb; and a 1 Mey lower 7~ level could only beexplained by an extremely 
large configurational interaction in this case. In addition there should be 
a 5~ level in the neighbourhood of the 9~ level. This 5~ level is not found 
in **Pb and only slightly below the 9~ level in Pb. Considering also 
the position of the 5~ level in ?°6Pb (see fig. 1) this suggests that for lower 
mass numbers the 5~ level might be still lower. In this case the half-life 
of the isomeric level would be determined by a rather fast E4 (9-  5-) 
transition, which may explain the fact that we were not able to observe 
isomers in 2°Pb and 1°8Pb. 


* This work gives references to previous investigations of 204™Phb. 
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Tl. Odd-odd isotopes —t}°*T1 has been shown (Bergstrém et al. 1953, 


_ Passel et al. 1954) to decay with three y-rays in cascade having energies 


261 (M4), 282 (M1+-E2) and 48-7 (E2) kev. In the decay of the 3-5 h 2°mPb 
_ activity three y-rays of energies 460 (E3), 389 (M1) and 401 (E2) kev 
converted in Tl were found. In both 1°*™T] and 2&mT the sequence of 
the non-isomeric y-rays has not been established. If the ground states of 
1981 and 2°T' are 2-, the excited levels are 2-, 3~ or 4-, 5-, 9+ and 2-, 3- 
or 4-, 5-, 8t respectively. In addition in the decay of 1°°Pb and ?°Pb 
strong H2 y-rays are found of 173 and 148 kev respectively, probably 
originating in 0~ states decaying to the groundstate. 


Fig. | 
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Levels of some Pb and T! isotopes. 


Tl. Even-odd isotopes.—In the isotopes 2037], 2017], Tl and 19’Tl 

the first excited states (presumably d;,, states) are situated 279, 330, 367 
and 388 kev respectively above the 84/2 ground states. The second 
excited states (presumably d;,, states) are 678, 692, 722 and probably 
774kev respectively above the groundstates. It is probable that the 
isomeric h,,/. levels are situated above the d;/. states and therefore do xe 
decay by M4 transitions to the dss states ae by He ee ay e 
d;/2 states, as is also the case in Au (Gillon e¢ al. 1954). The half-lives 
would then be short which may explain why only 197™T1 has been observed. 

al. 1953). 

ae eae . noticed Ane in the decay of 2°4™Pb envee y-rays are eusenad. 
According to the shell model only one or two should be ee WN 
cannot be excluded that one or two of the y-rays reported in * are 


emitted from 29!™TI. 
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IX. Neutron Deficient Isotopes of Pb and TI~II: Mass Numbers 
204, 202, 201, and 200 


By K. E. Brrexvist, I. Brrasrrom, C. J. HERRLANDER, S. HuLTBerG, 
H. Siatis, E. Soxontowsx1, A. H. Warsrra* and T. WIEDLINGT 


Nobel Institute of Physics, Stockholm, Swedent 
[Received September 20, 1954] 


THE investigation of the decay of 2°°-22Pb was started mainly in order to 
obtain more detailed information about ®°°"Pb (Maeder and Wapstra 
1954, Maeder, Wapstra, Nijgh and Ornstein 1954) and, since it has been 
suggested that 1°°Pb should also have a high spin isomeric state (Passel, 
Michel and Bergstrém 1954), to check whether a similar case of an even— 
even isomer existed in 2°°Pb. 

Tl was bombarded with deuterons (H,,,,=25 Mev) in the Stockholm 
cyclotron and with protons (Hy,,=190 Mev) in the Uppsala synchro- 
cyclotron. The almost carrier free lead activities were evaporated in 
vacuum on to thin Al foils (0-15 mg/cm?), which were used as sources in a 
two-directional focusing B-spectrometer (r=50 cm, resolution ~0-3°%). 
Conversion electrons were investigated in the energy region 10-1600 kev. 
The preliminary results of the measurements are given below. (A 
detailed report of the work still in progress will be given in Arkiv for 
Fysik.) The energy values given in table 1 have an accuracy of 0-2% 
or better. The K/L ratios may in some cases be uncertain. In spite 
of the high resolution used, the large number of conversion lines makes 
the overlapping of lines rather probable. 

204™Pb —Krohn and Raboy (Krohn and Raboy 1954) recently found by 
delayed coincidence measurements (Nal counters) that, in addition to the 
isomeric 913 kev y-ray in ?°*"Pb there existed another y-ray with a 15 kev 
lower energy. Figure 1 shows the K-conversion lines of the 913-0 and — 
899-3 kev y-rays. If the 913, 374 and 899-3 kev y-rays are emitted in 
cascade, as suggested by the experiments of Krohn and Raboy, the 
intensity ratio K go9.3/Ky,3=0'12 (expected value=0-12) suggests that the 
899-3 kev is due to a pure E2 transition assuming that the 913 kev y-ray 
is due to an ES transition. 

202™Ph gnd 202”T].—Due to the high resolution of the f-spectrometer, 
some uncertainties in the decay scheme suggested by Maeder and Wapstra 
(1954) could be eliminated. Some weak conversion lines assigned to 
possible y-rays from levels in °°’Pb were now shown to pe converted in 
Tl, suggesting a 10% branching from the isomeric state in ?%Pb to a 
ON es a ee ne 


* On leave from the Instituut voor Kernphysisch Onderzoek, Amsterdam. 
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Table 1 
Ke 
Nuclide EH, | conversion |Converted| Multi- Based on 
intensity in polarity 
201mPh 
68 min | 913-0 100 Pb E5 Earlier assignment 
899-3 12-:0+1 Pb E2 K899-3/K913=0-12 
202mPbh | 129-4 2 Pb K4 712; K/L<0-003 
3-5 h Ly/Lyy=1-7 
421-9 35 Pb K2 K/L=2°5 ; 
aK =0-035 
657-6 3°5 Pb El Kil i-5% 
. aK =0-005 
786-9 100 Pb E5 Tia, K/L=1-35; 
aK =0-089 
961-2 12 Pb K2 K/L=5:9 ; 
aK =0-006 
389-2 27 Tl M1+E2 ax =0-16 
401-3 5 Tl K2 aK =0-030 
459-9 14 Tl E3 K/L=1-51+0°15 


—_— 
bo 
eo) 
—_ 

| 


0-7 Only K-linje observed 
310-0 0:5 Tl EK2 ? K/L=0-7 2? 
330-3 100-0 nt M1+ E2 ? K/L=4-7 
361-2 25 Tl M1 K/L=5-6 
394:3 0-6 Tl K2 ? K/L=2:-5 
405-6 2-7 Tl E2 ? K/L=2:-8 
585 1-9 Tl K/L=5-5 
692 0-6 Tl E2 ? K/L=6-5 
708 2 0:3 Only K-linje observed 
766 1-0 Ti K/L=5-5 
825 ? 0-4 Only K-linje observed 
907 1-3 fi! K/L=5:5 
946 1-6 Tl K/L=6-5 
1099 0:2 A K/L=5:5 
30-6 Hg M1 L,/Ly,=10 
32-1 Hg M1 rots 
135-0 60 Hg Ml L,/Ly,=10 ; 
K/L=7-0 
167-2 100 Hg M1 L,/Ly,=10 ; 
K/L=6-2 


See 
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K- 
Nuclide Ey, conversion |Converted} Multi- Based on 
intensity in polarity 
200Ph 
21:5+0-4h]| 32:8 Only wieak M and | N lines observed 
109-5 3 alt M1 No L 
142-2 43 Tl M1 Kar6- 
no L 
148-0 100 Tl E2 K/L=0-4 ; 
Ly Ly 2350) 
: TIT 
158-9 4 L-line hjidden 
235-3 22 Tl M1 K/L=5-5 ; no Ly 
257°3 Cf Tl M1 K/L=4-5; no Lyy 
268-0 14 lil MI K/L=4:5; no Ly; 
289-5 il Wieak, no L-ljine observed 
450 Onlly M and N |lines observed 
2007] 
27h 116-3 5 Hg 
252-1 4 Hg 
289-1? 2 Only K-linle observed 
367-6 100 Hg K2 K/L=2-0 
578°8 7 Hg K2 ? K/L=4 
628-42 1 Only K-linje observed 
660-2? 1 Only K-linje observed 
788-4 1 Hg K/L=5 
828 8 Hg K/L=4 
1205 9 Hg K/L=4:3 
1224 1 Hg 
1363 1 Hg 
1515 1 Hg 


high spin isomeric state in 2°?Tl (odd-odd) from which three y-rays of 
energies 460, 389, and 401 kev (3, M1, and K2) may be emitted in cascade. 

201Pb and 2°"T1.—As can be seen in table 1, the decay of ?'Pb is very 
complex. In addition to the y-rays listed, there are several weak con- 
"version lines with a half-life ~ 9 hours, which are not assignable. By 
e- e--coincidence measurements in a double lens spectrometer* (Gerholm, 
to be published) the strong 330 and 361 kev y-rays were shown to coincide. 
Because of their intensities (cf. table 1) the latter y-ray should be emitted 
first. The 692 kev y-ray may be due to a cross over transition and the 
361 and 330 kev y-rays emitted from the first two excited states in ?"TI, 
which probably are d;. and d;. states. There were no conversion lines 
indicating a long-lived isomeric state in ave. 

It has previously been shown that the 167 and 135 kev y-rays from the 
excited levels in 2°'Hg are M1 transitions (Bergstrom, Hill and de Pasquali 


ES Se ae eT 
* We are indebted to fil. lic. T. R. Gerholm for putting his spectrometer at 
our disposal and for help with the measurements. 
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1953). The L,/Ly, ratios of the 30:6 and 32-1 kev y-rays are both ~10, 
showing that also these two y-rays are due to M1 transitions. e7 
y-coincidences in a f-spectrometer showed that the 30-6 and 32-1 kev 
y-rays both coincide with the 135 kev y-ray, but not with the 167 kev 
y-ray. No coincidences were found between the 135 and 167 kev y-rays. 
Assuming an sj. ground state in 2°!T] the following level assignment for 
201H@ agrees with the experiments: pj. (0 kev), 5/2-(1:5 kev), 3/2- 
(32:1 kev), and 1/2-(167-2 kev). No transitions were detected between 
the 167-2 and 1-5 kev levels, which may be explained by the low emission 
and conversion probability of E2-rays. Since all the observed y-rays 
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are M1 transitions the level scheme requires that the intensity ratio of 
the 30-6 and 32-1 kev y-rays is (30-6/32-1)?=0-86 and that of the 167 
and 135 kev y-rays (167/135)8=1-9. This is in agreement with the 
experimental values 1-0 and 3-2 supporting the suggested level scheme. 

2Pb and »°Tl.—e-—e~ coincidences in a double lens spectrometer 
(Gerholm, to be published) showed that the following y-rays coincide with 
the strong 148 kev y-ray in 2°Pb: 109, 142, 235, and 268 kev. The 
159, 257, and 289 kev y-rays do not coincide with the 148 kev y-ray. 
All y-rays but two (159 and 450 kev) fit into a level scheme in agreement 
with the measurements. The excited levels have the energies: 148, 
257, 289-5, and 525 kev. The parity and spin assignments are still 


uncertain. There was no evidence for an isomeric state in 2°Pb, having 
a half-life larger than about one hour. 
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The measurements on 2°T] are yet incomplete. However, the y-rays 
found confirm earlier measurements (Bergstrém, Hill and de Pasquali 
1953). There are several other conversion lines having a 27 h half-life. 
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THE present investigation of Pb and TI isotopes with A<200 is closely 
connected to similar work on higher mass numbers going on at the Nobel 
Institute in Stockholm (cf. I and II). Mainly because of the short half- 
lives in this region, however, the measurements have been performed at 
the Gustaf Werner Institute for Nuclear Chemistry. 

The activities, produced by TI] (p, xn) reactions in the Uppsala synchro- 
cyclotron, have been studied in a two-directional focusing 6-spectrometer 
(r=18-5 cm, resolution ~0-3%) (Arbman and Svartholm). Most spectra 
were obtained by means of an automatic device, using a ratemeter 
and a Speedomax recorder. For mass assignments and in some cases 
when preparing f-spectrometer samples, an electromagnetic isotope 
separator of the 90° type (Andersson ef al.) has been used. 

So far the conversion electron spectrum has been investigated only 
in the region 50-700 kev. All the activities except "Pb were shown 
to decay by electron capture. A full account of the investigations, still 
in progress, will later be published in Arkiv for Fysik. The following is 
a summary of the preliminary results (cf. table 1). 

!=™Pbh and Pb.—A 12 minute activity has been assigned to A=199 
partly because of excitation relations. Moreover, the decay curve of a 
mass-separated sample of A= 199 showed a strong ~10 minute activity. In 
the same way a 90 minute activity was assigned to the same mass number. 
Since this activity was not found in the electron spectrum of a pure Tl 
sample and the conversion definitely occurs in Tl, the 90 minute activity 
must be due to the ground state. It is therefore concluded that 1°Pb 
has an isomeric state, not known earlier, agreeing with the fact that the 
electron lines of the 12 minute activity are converted in Pb. 

If the 354 and 367 kev y-rays are emitted in cascade from levels in 
Tl as indicated by a possible 722 kev y-ray, the 367 kev y-ray should 
be emitted from a lower level for intensity reasons. 

“"Pb.—There are indications that the 169 and 234 kev y-rays are 
converted in Pb, suggesting an isomeric state which may have about the 
same half-life as the ground state. (The half-lives of the y-rays in the 
table have not all been measured with the given accuracy.) 


* On leave from the Instituut voor Kernphysisch Onderzoek, Amsterdam. 
+ Communicated by Professor Siegbahn and Professor Svedberg. 
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Table 1 


Noorae E, K-conversion| Converted | Multi- 


(Ly + Ly)/Lyyr=1 4 
MI K/L2:5, no Ly 
M1 K/L&5, no Ly 


M1 K/LA6 


196] 2 


~4 h 


ee in polarity ESO 
199mPh 
12-2+0-3m | 423 Pb M4 ? K/L=1-9, no 
y-ray detected 
199Ph 
90+10 m M1 K/L&5, no Ly; 
M1 K/L5, no Lyyy 
198Ph 
2:3+0-2h K2 K/L=0-5 
(Ly +-Lyp)/Lr=1-9 
M1 ? No Ly; 
2 K/L0-4 
(Ly+-Lyy)/Lyyy=2°3 
K2 ? K/L=0°3 
M4? | (Ly +-Lyy)/Lmr=2-0 
M1 K/L&5, no Lyyy 
MI K/L=5, no Ly; 
1977] =" 
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199T] and 1°8T]—The y-rays found in the decay of these nuclides 
confirm earlier measurements (Bergstrém et al. 1953). As to 19Tl 
the y-rays reported by Bergstrém et al. (1953) and Passel e¢ al. (1954) 
were not observed. One possible explanation is that the isomeric state 
of 198Pb has not the half-life ~25 minutes as suggested by Passel et al. 
(1954), but is so short-lived that it decays considerably before the chemical 
separation of Pb from the TI target. 

197T|_—The 134 kev y-ray found in the decay of 1°’Tl is undoubtedly 
identical with the y-ray of the same energy emitted in the decay of 
re: 

196T] —The assignment of the 426 kev y-ray to 1°°Tl is supported 
by the following facts. Excitation relations indicated A< 198. In the 
electron spectrum of A197 (mass separated), however, the conversion 
lines of the 426 kev y-ray did not show up, in spite of their strong 
intensities. Furthermore a 426 kev y-ray has been observed in the 
B-decay of 19®Au (Staehelin 1952). 
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Cubic Metals 


By A. J. E. Foreman and W. M. Lomur 
Atomic Energy Research Establishment, Harwell, Berks.* 


[Received October 2, 1954] 


In interpreting the modes of deformation of a metal crystal in terms of the 
theory of dislocations it is usually assumed that the direction and plane of 
slip are determined by the criterion of minimum dislocation energy. 
Since the energy of a dislocation varies as the square of the Burgers 
vector, it follows that the slip direction is the closest packed direction in 
the crystal, in agreement with observation, but a more detailed analysis 
appears to be necessary. before the occurrence of slip on only certain 
prominent crystallographic planes may be adequately explained. 

Calculations have recently been made of some dislocation energies in the 
face-centred cubic lattice, using the anisotropic elasticity theory of 
Eshelby, Read and Shockley (1953), for various edge, mixed and screw 
dislocations in several different slip planes containing the (110) Burgers 
vector. Eshelby (1949) has shown that the energy per unit length of a 
straight dislocation line (Burgers vector b) takes the form 


Kee, R 
1) b= Se aa teal omer) 


—— nh 
2 hie Us 
where K is a function of the elastic constants, edge-screw composition 
and dislocation slip plane ; & is the radius of the dislocation strain field 
(i.e. the mean spacing in a dislocation array, typically about 10*b), 
and ry is an effective radius of the central core of the dislocation. For 
elastic isotropy the energy factor K becomes independent of the slip 
plane orientation and is given by 
sin? 6 
K=6 (ost 0+ =): Pe ee Be OF 
where @ is the angle which the Burgers vector makes with the dislocation 
line, G is the shear modulus and v is Poisson’s ratio. 
In the numerical calculation of K we use the following values for the 
single crystal elastic constants (units of 10” dynes/cm?). 


Oe Ce On A Reference 
Al 1.08 0-622 «0-284. 12 Schmid and Boas (1950) 
Au 187 1-57 es 


Cu 1-710 1-239 0-756 3:2 Lazarus (1949) 


The ‘anisotropy factor’ A=2C44/(Cyi—C12) is unity for an isotropic 
crystal. 


* Communicated by the Author. 
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The calculated values of the energy factor K for aluminium and copper 
are shown in the figure for the {100}, {110} and {111} dislocation slip 
planes, with curves of best fit drawn through the calculated points. The 
curves for gold fall near those for copper and have been omitted from the 
figure for the sake of clarity. The notable features of these results are 

(a) the variation of K with @ still approximates to the isotropic form 

K=K (screw) cos? 0+ K (edge) sin? @ leet 
for each slip plane, 

(b) except for a pure screw dislocation the spread of values for different 
slip planes increases with increasing anisotropy, 

(c) the ratio K (screw)/K (edge) has the order of magnitude of 1—», 
where v is the mean (polycrystalline) value of Poisson’s ratio, 

(d) for none of the metals considered has the observed {111} slip plane 
corresponded to the lowest K values. 


ES 
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Edge-screw composition angle 0 
beers of A with the angle @ between the dislocation line and 
(110) Burgers vector, for three dislocation slip planes. 


‘The value of 7) in eqn. (1), also, will depend on the choice of dislocation 
slip plane, especially since dissociation into a pair of partial dislocations 
connected by a stacking fault may occur in the {111} plane, with a con- 
sequent lowering of energy. Seeger and Schdck (1953) have investigated 
this extended dislocation, using the Peierls dislocation model with aniso- 
tropic elasticity theory and taking as the stacking fault energy the coherent 
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twin boundary energy of the crystal. For copper the stacking fault 
width is shown to be about 8 b for a screw and 12 b for an edge dislocation, 
and the corresponding energy gains as compared with the energy of a 
narrow Peierls—-Nabarro dislocation may be shown to be about 8% and 
10% respectively, if we take the dislocation strain field radius as of the 
order of 10*b. In the case of aluminium the stacking fault is only 1 or 2 
atoms wide, which hardly constitutes a true dissociation, and the resulting 
energy gains are only a few per cent. 

These results indicate that even when the gain in energy due to dissocia- 
tion is taken into account the energy of an edge dislocation in a {111} 
plane still exceeds that for a {110} plane. The screw dislocation is, how- 
ever, energetically favoured in the close packed {111} plane by a factor of 
about 8°% for copper and a few per cent for aluminium. Thus the mean 
energy of a dislocation ring in a {110} plane is approximately equal to that 
for a ring in a {111} plane, and it therefore seems unlikely from energy 
considerations that the dislocation content of a crystal as grown should be 
composed of predominantly {111} dislocation rings. 

It is suggested that dislocations may initially exist on numerous planes 
in a face-centred cubic crystal, but that the pure screw segments of any 
general dislocation loop will become extended into a close packed {111} 
plane and their subsequent glide motion may then be limited to the plane 
of dissociation. It is not at present clear what is the smallest degree of 
dissociation which is sufficient to effectively capture the screw dislocations 
into the {111} plane. Ifthe remaining mixed and edge-type segments are 
sufficiently mobile to glide under the applied stress they could give some 
slip on planes other than {111}, but it is unlikely that this would be in 
observable amounts unless some dislocation multiplication mechanisms 
were to operate. A Frank—Read source is not expected to work on any 
plane other than {111}, however, since it would require the expanding loop 
of dislocation to repeatedly pass through the screw orientation without 

i tured into a {111} plane. | 
ee factor ttn vee aan to the choice of slip plane is the 
lattice mobility of the dislocations. Nabarro (1952) has suggested that 
when a dislocation becomes extended the decrease in elastic energy is 
accompanied by an increase in mobility, so that dislocations would be 
expected to glide more easily on the {111 } planes of a face-centred cubic 
metal. Aluminium is an interesting metal in this respect, since the nearly 
isotropic elastic properties and small dissociation give the dislocation 
energy only a small dependence on the slip plane orientation. The 
dislocations in the {111} plane will, however, be slightly wider so that 
dislocation mobility, which is a sensitive function of width, is ee pooted to 
be greatest in this plane. This could explain the BOER. of {111} 
slip at room temperatures ; the occurrence of slip on. oes. pee at 
higher temperatures is in agreement with the conclusion that the depen- 
dence of dislocation energy and mobility on slip plane orientation is 


; A J cubic metals 
smaller for aluminium than for most other face-centred cubic metals. 
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XII. Soft X-Ray Spectra of Magnesium—Aluminium, Magnesium- 
Silicon and Aluminium—Silicon Alloys 
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ABSTRACT 

The soft x-ray L spectra of the elements Mg, Al and Si have been 
determined in all three of the sets of binary alloys, and have been com- 
pared with the spectra of the pure elements. Several differences appear : 
(a) a widening of the emission edges with alloying, (b) in one case a 
probable diminution of band width and (c) structural changes in the shapes 
of the bands. We can attempt to understand these effects on the 
assumption that the alloys consist of two or more phases. No case of 
identity of the bands of the two constituents, such as might be expected 
from a perfectly ordered single-phase alloy, has been found. 


$1. InrRopucTiIon 


THE L, 5 emission spectra of magnesium, aluminium and silicon in 
binary alloys of magnesium—aluminium, magnesium-silicon and alumi- 
nium-silicon of different percentage compositions have been studied 
with a ruled grating grazing incidence spectrograph. An investigation 
of the K spectra of these metals in some alloys has been made by Das 
Gupta (1946) and in magnesium—aluminium alloys by Farineau (1938). 
The L, 3 spectra of magnesium, aluminium and silicon in a pure state 
have been investigated by Skinner (1940), and the structures in the 
bands have been interpreted in terms of the structure of the electron 
bands in the solid state. The high energy kinks in the spectra are 
evidence of the overlapping of Brillouin zones in these metals. Our 
purpose is to find out how the structures of the bands are changed on 
alloying. 
§2. EXPERIMENTAL TECHNIQUE 


The spectrograph was the same as that employed by Skinner (1940). 
The concave glass grating of one metre radius of curvature having 
30 000 lines per inch was set at 6° grazing incidence and the region 
eee ee arma Es 

* Now at the Physics Department, University College of Science, Calcutta 

+ Communicated by H. W. B. Skinner, F.R.S. 
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40-400 a could be investigated in one fixed setting of the spectrograph. 
The slit system consists of three slits S,, S, and 8, set in proper alignment 
so that about 15 mm of the ruled grating are utilized to produce the 
spectrum. The slit S, immediately in front of the anti-cathode which 
acts as the source was 0-03 mm in width. The middle slit S, is about 
1-5 mm and the slit S, nearest to the grating is about 0-7 mm. 

In the designing of the x-ray tube, the technique of Skinner (1940) 
has been essentially followed. The modification that has been made 
consists mainly in introducing a liquid air trap very close to the copper 
target. A glass bulb of one litre volume was used for the x-ray tube. 

It is considered preferable to use a glass tube rather than a metal tube 
from the point of view of cleaning and degassing, and also for observation 
during the exposure. The anticathode was a square sectioned copper 
tube, of 8 mm side and } mm wall, mounted vertically on a glass cone. 
In this way any of the 4 faces could be used as a target for the exposure. 
If radiation is not to be wasted, the position of the anticathode face 
must be adjusted so as to present an oblique face to the slit S, which lies 
within the acceptance angle of the spectrograph. This involves the 
accurate adjustment of the axis of the anticathode relative to the spectro- 
graph slit. In order to do this, a flange on the glass bulb is sealed by 
means of wax to the flange on which the slit is mounted, so that, by warm- 
ing the wax, the bulb may be slid relative to the slit. The glass bulb is 
joined to the pumping system by a copper bellows in order to provide 
latitude for this adjustment. The electron gun was a vertical coil of 
nickel wire with oxide coating. This was placed a few mm from the 
anticathode face, and was surrounded by a cap in order to concentrate 
the beam on to an area of about 2x10 mm on the anticathode. An 
electron current of 50-100 ma at 3 kv could be obtained. This method, 
using a dull emission filament, rather than the tungsten filament used by 
Skinner, is convenient, but has the disadvantage that it is very difficult 
to avoid all traces of a Ba line at 165 a which may in certain cases be a 
nuisance. 

It is most important to ensure that the surface of the metal used as 
target is free from contamination (oxide, ete.). For pure metals, there- 
fore, we adopted the method of Skinner (1940) and evaporated thin 
layers of the material in a high vacuum on to the copper anticathode. 
For this purpose, an evaporating unit was placed behind the anticathode, 
on the opposite side to the electron gun. The problem for alloys is more 
difficult, and three methods were used. 

(1) The two metals A and B to form the alloy are placed in small 
lumps separately in two conical tungsten filaments, in front of the 
target of the x-ray tube. When the vacuum of the order of 10-6 mm 
is attained, the filaments are raised to two different temperatures by 
regulating the current passing through them to adjust the rates of 
evaporation of the metals A and B so that an alloy of suitable com- 
position is deposited on the target by the simultaneous evaporation of 
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the metals A and B. The evaporation rate of the metals must be slow 
and steady in order to give a homogeneous deposit. When a sufficiently 
thick layer has been obtained, the deposition must be ended by turning 
the target with the evaporation still running in order to ensure that the 
surface layer has the right composition. The target with the thin 
deposit of the alloy was then heated for about one hour at about 400°c 
in an attempt to make the structure homogeneous and ordered. It is 
difficult to judge how effective this procedure was. This method has the 
disadvantages that the composition of the alloy cannot be pre-determined 
with any accuracy, and also that the percentage can only be very roughly 
determined from the exposure time and density of the spectra obtained. 

(2) The alloy is prepared beforehand in an induction furnace and the 
percentage composition is accurately determined by chemical analysis. 
The prepared alloy is placed in small lumps in a conical tungsten filament. 
When the proper vacuum is attained, the tungsten filament is raised very 
quickly to a very high temperature so that the entire lump is immediately 
evaporated and a deposit on the water-cooled copper target is formed 
having the same composition as that of the starting material. Several 
photographs were taken of one particular composition and the repro- 
ducible results justify the assumption that the composition was approxi- 
mately maintained. 

(3) A rectangular strip of alloy 2 cm x5 mm of thickness varying from 
0-3-0-5 mm of known composition was fixed to the water-cooled rect- 
angular face of the copper target. A good contact of the alloy strip 
with the flat copper face of the target was obtained by grinding the alloy 
strip, as well as the anticathode, flat. The strip of alloy clamped to the 
target: was cleaned in vacuum at regular intervals during the exposure 
by means of a steel wool brush. This method of obtaining a surface 
free from impurity turned out to be surprisingly effective, and of course 
the constitution of the alloy is much better guaranteed than with methods 
(1) and (2). But the method has the disadvantage that, since there is 
only a poor heat contact between the specimen and the copper anode, 
high excitation currents cannot be used. 

- Some of the alloys used in methods (2) and (3) were kindly supplied 
by the British Aluminium Company, with a report of the chemical 
analysis. 

For methods (1) and (2) the x-ray tube was usually run at 3000 volts, 
50-75 ma, and for method (3) 3000 volts, 20-30 ms. Exposures ranged 
from about } hour for the pure metals up to several hours for weak alloys. 

Several standard photographs of L, edges of pure magnesium and 
aluminium were taken in Ist and 2nd orders. The sharpness of the 
emission edges (of the order of 0-20 ev at room temperatures) was used 
as a check on the proper focusing of the spectrograph. Photographs 
were taken with different periods of exposure and from the photometer 
curves the most suitable exposure period to give a good photograph of 
a band was estimated. 
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Ilford Q, emulsions coated on to thin glass about $ mm thick were 
used throughout. No method of calibrating the blackening of the plate 
with intensity was used, and this must be borne in mind in studying the 
results to be given, because the denser exposures show definite signs of 
saturation. 

The dispersion curve was obtained from the known wavelengths of a 
characteristic point on the K-band of graphite, which, by painting 
colloidal graphite on to a face of the anticathode, could be obtained in 
8 different orders. These were found to agree well in form and wavelength 
with those obtained previously by Skinner (1940). Three spectra could 
be taken side by side by the use of a shutter in front of the plate, and it 
was usually arranged that the spectrum of the pure metal was taken side 
by side with the spectrum of the element in an alloy being investigated. 
In fact, if we are using method (1) of alloy preparation, the juxta-position 
of a spectrum of the pure metal is necessary so that one can judge roughly 
the composition of the alloy by the relative intensity of the two bands. 


§ 3. EXPERIMENTAL RESULTS 
3.1. Mg—Al Alloys 


In fig. 1, the L, , spectra of Mg and Al are shown for various com- 
positions of alloy. The corresponding spectra for the pure metals are 
also given for comparison. 


(a) Sharpness of Emission Edges 


Tn all the curves there is a sharp drop at the short wavelength limit 
of the band. This feature is known as ‘the emission edge’. It is seen 
clearly separated into two components, the L, and L; edges, correspond- 
ing to the two x-ray levels excited. In all cases, there is no appreciable 
shift in the energy of the edge due to alloying. The L, edge breadth 
may be defined as the energy gap between the beginning of the drop in 
intensity of the L, edge and the beginning of the corresponding drop in 
the L, edge. Table 1 shows the results. The alloy compositions which 
are asterisked were determined chemically and were used as strips fixed 
to the target (method 3). The remainder were obtained by simultaneous 


evaporation (method 1), and therefore the compositions are only roughly 
known. 


Table 1. Observed L, Edge Breadths for Mg and Al in Mg—Al Alloys 


HAL Al edge breadth Mg edge breadth 

ev ev 

; Oe Cue 

5 _- 0-4 

9 a 0-4 

50 0-7 0-3 

60 0-5 0-4 

hes 0-8 0-5 

93* 0-6 a 

95* 0-7 — 
100 0-2 


: 
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The values for the pure metals agree with those obtained at room 
temperature by Skinner (1940). 


It will be seen that there is a marked rise in the edge breadths both 
for Al and Mg with even an alloy containing 5% of the other metal. 


Fig. 1 
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A higher percentage does not appear, within the limits of error, to have 
much further effect. he . 
The increase is a very clear effect. In a pure metal, Skinner (1940) has 
shown that the edge-widths correspond fairly accurately with the statisti- 
cally calculated breadth of the Fermi edge at different temperatures. 


G2 
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But in the alloys, the edges are considerably wider than could be 
explained as a temperature effect ; in fact the alloys would have to 
be molten to give so large a breadth. One may perhaps try to think of 
this effect as due to slight differences in the energy of individual transitions 
for emission by atoms differently placed in the lattice. In a pure metal, 
of course, each atom is surrounded by an identical array of atoms, but in 
an alloy, different arrangements of neighbours are obviously possible. 

It is to be observed that we did not get any case of complete disappear- 
ance of the edge, such as was reported by Yoshida (1935) for the Al K 
band from a dilute solution of Al in Cu and by Skinner and Johnston 
(1937) for the Al L band and Be K band for weak solutions in Cu, and 
also, in the latter case, in Al. But it must be noted that we have not 
in the present experiments on Mg—Al alloys been able to use an alloy of 
less than 30°, whereas the effect reported by the other authors only 
appears in the Al-Cu alloys at less than 10%, and in the case of the Be 
alloys at roughly 1 to 2%. It is possible that the disappearance of the 
edge only occurs in the case of genuine solid solution. This cannot be 
proved from Skinner and Johnston’s results on account of the considerable 
inaccuracy of their estimate of composition, but there certainly seem to 
be signs of the development of edges as the percentage rises to a figure 
much beyond the solubility limits. For comparison, these are given in 
table 2. 


Table 2. Approximate Solubility Limits (% by weight) at Room 
Temperature for Various Alloys 
Alin Mg Alin Cu Mg in Al Be in Cu Be in Al 
2 9 2 0-2 2 
(b) Band Shapes 

The shapes of the bands for pure Mg and Al also correspond well with 
those of Skinner (1940). In Al, there are two clear kinks at B and C, 
but in Mg only one. 

Figure 1 shows that the kink C in Al is obliterated on alloying. It 
is, however, possible that this is simply related to the widening of the 
edge. A similar obliteration was shown by Skinner to occur on simply 
heating pure Al to 670°K, in which case the edge breadth was 0-4 ev. 

There is no clear movement of the point B on alloying in the case of 
Al, but there is certainly a definite progressive movement towards lower 
energies in the case of Mg. We might interpret this as an increase in 
the zone overlap due to the increased number of conduction electrons 
in the band on account of the fact that Al has 3 valence electrons as 
against 2 for Mg. Such a simple argument would lead to a movement of 
B for Al to higher energies on alloying. Actually B stays nearly fixed, 
and so the argument is probably oversimplified. 


(c) Band Widths 


With regard to the band widths, although many of the traces of fig. 1 
are not suitable for accurate estimation, it will be seen that the L-bands 
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for both Al and Mg are not greatly altered in width by alloying. The 
Mg band is about 7 ev in width and the Al band about 11 ev, as com- 
pared with Skinner’s results for the pure metals, namely 7-2-+0-2 and 
13-2+0-4. Our band width for pure Al therefore does not exactly agree 
_ with Skinner’s value, but the low energy end of the photometer trace is 
not satisfactory for an accurate estimate, and subsequent work by Jacob 
and Wood has confirmed that Skinner’s width is correct. There seems 
therefore to be a small decrease in the Al band-width on alloying with 
Mg. This is what would be expected if we assume that the band corres- 
ponds to an alloy with an average of rather less than 3 electrons per 
atom. The Mg bands are not sufficiently good to judge whether there 
is a corresponding increase in the band width of Mg. It is clear that 
there is no equality of band widths of Al and Mg in any alloy. Farineau 
(1938) reported such an equality in the case of the K-bands of Mg and Al 
in the compounds Mg,Al, and Mg,Al,, which gave band widths of about 
10-5 volts in each case as against his results of about 8 and 13 volts for 
the pure metals. The band shapes also changed into an identical form 
for the two constituents, possibly slightly different for the alloys Al,;Mg, 
and Al,Mg;. This result is certainly what one would expect in the case 
of an alloy which is a perfectly ordered inter-metallic compound. Fari- 
neau remarks that it is difficult to obtain such an alloy in a homogeneous 
phase. So, although we have covered the range between these two 
alloys (37:5% to 57-5 Mg by weight), and although we used heating to 
400°c on the target for ordering, there is no positive test that our alloys 
were homogeneous. In view of the importance of the matter, further 
work is required, with an alloy of proved homogeneous phase. 

At present it is clear that, at any rate for alloys which are not perfectly 
ordered intermetallic compounds, the electron band in the neighbourhood 
of one constituent is not of the same shape as in the neighbourhood of 
the other. From the work of Skinner and Johnston (1937), this con- 
clusion also follows for dilute alloys of Be in Al and Cu. Since the 
band widths of pure Be, Al and Cu are not very different, no conclusion 
can be made about the equality of band widths in these cases. It appears 
from the present work for a whole range of Mg—Al alloys that neither 
the band widths nor shapes are the same, and, as we shall see subse- 


quently, the same applies to other alloys. 


3.2. Al-Si Alloys 

Our results are shown in fig. 2, which consists of tracings of photo- 
metric curves. The dotted parts of the curves are uncertain. The 
alloys were of known composition and were clamped to the anti-cathodes 
as strips, the surface being cleaned by the wire brush. Other eerie 
were obtained using the method (2) of evaporation, but ee similar 
results. The broadening of the Al L, edge is somewhat less pronounced 
than in the case of the Mg—Al alloys but as the edge breadth in the ne 
of the 5% alloy is about 0-5 ev, the effect seems present. It will be note 
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that the kink C is still clear in spite of the edge widening. As in the case 
of the Mg alloys, the kink B remains, but a new marked kink at A has 
appeared, which might be interpreted as additional zone-overlapping 


through alloying with tetra-valent silicon. 
Fig. 2 
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L, 3 spectra from Al-Si alloys. 


The Si band also appears to show structural changes in the case of the 
5% alloy. The band widths of both Al and Si seem little affected by the 
alloying, even though the percentages of Si are low. Our curve is not 

suitable for finding the exact value of the Si band width, but Skinner 
(1940) gives 18-2 volts, and the Al bands are definitely considerably 
narrower than this. The solubility limit of Al in Si is 24% and there is 
practically no solubility of Si in Al. 


3.3. Mg—Si Alloys 


Three alloys were investigated which were prepared by the method 
(1) of simultaneous evaporation, followed by heat treatment on the 
target. The proportions of Mg and Si could roughly be estimated from 
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the relative intensities of the respective L bands. The results are shown 
in fig. 3. The 50% alloy may perhaps show an additional kink in the 
Mg band at C, though this is not certain. The weak alloy gave only 
a faint trace, and the results are perhaps not to be relied upon. But as 
far as the evidence goes, a considerable extension of the band towards 
low energies is indicated. It may be that there we have an example of 
the equalization of the band-widths of the two component metals in the 
case of a weak alloy. However, the solubility limit of Mg in Si is very 
small. 

The Si band, in the alloy of 50% Si seems to show a definite structural 
change from that obtained with pure Si. 


Fig. 3 
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L, spectra from Mg-Si alloys. 


§ 4. CONCLUSIONS 


It seems clear that there is no one generally applicable picture of the 
effects of alloying on the soft x-ray bands. Three types of cases have 


to be considered :— pee | 
(1) There is some evidence by previous authors that in the case of 


genuine solid solutions, the edge of the band of the dilute pate 
disappears completely. In these cases, there are phate ee e ee na 
which take part in conductivity in the neighbourhood of the fe 4 ie 
dilute constituent. These atoms may perhaps act as onors’ 0 
electrons to the conduction band of the other atoms. 


86 On the Soft X-Ray Spectra of Mg—Al, Mg-Si and Al-Si Alloys 


(2) There is the evidence of Farineau that in the case of the inter- 
metallic compounds Mg,Al, and Mg,Al, the K-bands of both constituents 
are identical in shape and width. Although we have not worked with 
alloys of these definite compositions, in the case of the L-spectra we 
have not been able to find any signs of a similar case. Further work is 
needed to confirm Farineau’s conclusion. 

(3) There remain the rather nondescript set of ordinary alloys to which 
most of the present work applied. These consist of arbitrary proportions 
of the constituents and in these cases it may be that heat treatments, 
etc. are important factors. We have shown that in one case at least 
(Al-Mg alloy) the band-width of the Al is probably slightly diminished, 
though in the other cases the evidence is not good enough to draw any 
conclusion. We have also shown that the bands undergo structural 
changes on alloying, and that the edge-widths are broadened. These alloys 
must of course be expected to consist of mixtures of different phases 
and the results could be interpreted on the assumption that the alloys 
exist in the form of clusters of atoms of each type with a proportion of 
atoms of the other type included, rather than as a uniform mixture of 
both types of atom. 


The authors acknowledge gratefully the interest taken in this work 
by Professor H. W. B. Skinner, F.R.S., who suggested it and was also 
partly responsible for writing up the results. 
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SUMMARY 


A plastically deformed specimen of nickel is assumed to consist of an 
assembly of spherical ‘ grains’ each of which contains a uniform distri- 
bution of point defects. On this model, a calculation is made of the rate of 
evolution of energy when such a specimen is heated at a uniform rate. 
The result is compared with experimental data given by Clarebrough, 
Hargreaves and West and it is shown that the model adequately explains 
a part of the evolution of energy at low temperatures. The activation 
energy for migration of the defects is calculated to be 1-0 ev (23 kcal/mole). 
It is suggested that the defects are probably vacancies but the possibilities 
that they are interstitials or vacancy pairs are not excluded. On the 
vacancy hypothesis, the energy stored per vacancy is about 1-6 ev and the 
size of the ‘ grains ’ is of the same order as the size of the particles described 
by Gay and Kelly. 


§ 1. INTRODUCTION 


WuEN a plastically deformed specimen of nickel is heated at a uniform rate, 

a considerable fraction of the energy stored during the deformation is 

released before recrystallization. A typical curve for the rate of release 

of energy as a function of temperature, as obtained. by Clarebrough, 

Hargreaves and West (1953), is shown in fig. 1. There are two distinct 
Fig. 1 
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peaks, one at about 250°c and the other at about 600°C. These authors 
suggested that the release of energy represented by the peak at the lower 
temperature arises from the disappearance of vacancies created during the 
deformation. In the present paper this suggestion is put on a quantitative 
basis and it is shown that this peak can be explained as being due to the 
disappearance of an originally uniform array of point defects and that 
these defects are probably vacancies. It is further shown that the peak 
at the higher temperature cannot be explained on the same mathematically 
simple model. 
§ 2. THEORY 

As a basis for calculation, the following assumptions are made : 

(1) A deformed specimen of nickel consists of an assembly of spherical 
‘ grains ’ of radius a ; 

(2) after deformation, but before heating, each ‘grain’ contains a 
uniform distribution of point defects, their concentration being Cy ; 

(3) the defects move independently of each other by a diffusion- 
controlled process having energy of activation, Q, and frequency factor, 
Dy; and 

(4) the ‘ grain boundaries ’ act as infinite sinks for the defects and an 
energy, 1, is evolved whenever a defect reaches a boundary. 

If the whole volume, V, of the specimen is occupied by these ‘ grains ’, 
we have 3V/47a® ‘ grains’ and the rate of evolution of energy, 4P, is 
given by 

3V dN 
~ 4nat* * dt’ 
where N(t) is the total number of defects within a ‘ grain’ at time tf. 
Further, 


a= 


a 
N=4n | r2c(r, t) dr, 
0 
where r is a radial co-ordinate within any ‘ grain’ and c(r, t) is the 
concentration of defects at a point within a ‘ grain’. 
Therefore, we can write 


3HV (* oc 
= 1 ype 
oe as jr sh ae oso. tees A) 
and c(r,t) must satisfy the equation 
1 oc lo dc 
—_— = Vey So eS Sipe crs 
Di) tee pee (: *), eA et fh) 


where D(t)= Dy exp (—Q/RT(t)), T(t) being the absolute temperature. 
The solution of (2) subject to the appropriate initial and boundary 
conditions, viz. c(a, t)=0 for {> 0 and ¢(r, 0)=c, for 0<r<a, is 
2g 2% (—1)" . nar na? 
Se ea is >, 7 Sin —— exp | - ae D (7) dr | = eS) 
(cf. eqn. (95) on p. 29 of Barrer (1941) where the r in the denominator 
before the 2’ has been inadvertently omitted). 
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By substitution from (2) we can simplify (1) to the form 


3HV dc 
AP=— - D(t) E * |. pas, HED 
and on substituting (3) into (4) we obtain finally 
6Hz (oo) 2-2 ct 
Ap= A *% ni » exp | - =| D(r) ar |. eek (5) 
n=1 0 


Therefore, provided we know the dependence of 7' on ¢, we can in principle 
calculate 4P as a function of time. 


In the work of Clarebrough ef al. the temperature is raised at a uniform 
rate, «, so that we can write 
T(t)=T pat. (6) 


Using (6) and introducing a new variable u(t)=Q/RT, the integral in (5) 
can be written 
t D U —2% 
[ Pocxp [-Q/R(Tytorydr=— Mf" 
0 


aR J onr, © 


= 22 pu)—BQ/RT JL. - - » « (1) 


where f£(u)=Ei (—u)+exp (—u)/u. 
Then, defining 4P’=6HVc,D,/a2 and K=7?D,Q/a’«R, we obtain 
finally, 


AP= AP’ exp (—u) z exp [—n?K(B—fp)], . . - (8) 


where f, is the value of 8 for 7=7', i.e., at zero time. Thus, for given Q 
and K the problem of plotting 4P/4P’ as a function of ¢ can be reduced to 
the evaluation of 8 and then the summation of the series 


x exp (—n*y), 
ih 


where y=K (B—B,). 

The evaluation of B can be carried out by use of the tables and formulae 
given by Coulson and Duncanson (1942) or by use of the asymptotic 
series 


De aes} 
Bu) =Bi(—u)-+e-"fur=(e"p2)( 1 -- Fp oe ; aes elt) 


It is clear that f is a rapidly decreasing function of % 80 that, provided 7' 
is not very close to 7g, By can be neglected in comparison to B. nea ae 
In order to obtain 1% accuracy in the summation of the series, it is 
necessary to consider only the first term if y> 1-535 while two terms ae 
for y>0:57. However, for smaller values of y, corresponding to ci 
temperatures, the series converges only slowly and. it 1s convenient Oo 
transform it according to the relation (Whittaker and Watson 1946) 


co Teas ue =) a oe fn teaey. ¥ *2¥10 
Berra) 


y n=1 
n= 
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If only 1% accuracy is required the series on the right-hand side of (10) 
can be neglected in comparison with the first bracket for any y< 1-49. 
Thus it is simple to evaluate 1P/4P’ over any range of w. 

By differentiating the first term of (8), we find that this expression has 
one maximum, this occurring at a value of wu for which AB<1 ie. at a 
point where this term is not a good approximation to 4P/ AP’. However, 
on substituting the first bracket of (10) for the series in (8) and differ- 
entiating, we find that this expression also has one maximum value which 
occurs when 


K(B—Bo’=7 ((B—Bo)?—ke "wy. . ess CY 
By neglecting 8, and using the asymptotic series (9) for 8, (11) reduces to 
K= jure (: ¥ iia oe ee es 


From (12) and (9) it is clear that, at this maximum, KB~7/4<1-49. 
Therefore, 4P/4P’ should have a maximum at this value of wu and, from 
the discussion above, this should be the only maximum of 4P/4P’. 

Finally, by taking logarithms of both sides of (12) and substituting for 
K and u, we obtain 


47D,.R * 2 2 
jee ee =in| z 0-3 +.) _ (13) 


aQ  RT* T* Q Q? 
where 7'* is the temperature at which 4P is a maximum. 

Therefore, from a series of values of « and 7, i.e., of heating rate and 
temperature for maximum JAP, we can plot In («/7'**) against 1/7’* and 
hence obtain a first approximation to Q/R by calculating the slope by the 
method of least squares. The method of successive approximations can 
then be used to refine this value and to evaluate D,/a?, the linearity of the 
final relationship giving a check on the validity of the theory. When we 
have values of Q/R and D,/a? we can use the procedure outlined above to 
plot 4P/4P’ as a function of 7’ (or of ¢) and this curve can be compared 
with that determined experimentally. From the maximum value of 
AP.x >t, AP’ can be calculated and hence a value obtained for Eco. 

One relaxation of the basic assumptions can be made. Instead of 
assuming that all the ‘ grains’ have the same radius, a, we may suppose 
that their radii form a statistical distribution of fairly narrow width about 
some mean value. A recalculation on the basis of this model leaves the 


whole theory unaltered except that the a appearing in the definition of K 
must be interpreted as the ratio 


second moment (about zero) of the grain radii 
ae pt ee 
first 


29 ( 29 92 ) 9 2? 29 22 
while the a? in the definition of 4P’ must be interpreted as 


third moment (about zero) of the grain radii 
first ; : 


by ( 99 +P) ) 99 9 29 29 
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Therefore, the determinations of Q and K and of the shape of the final 
curve will remain unaltered. Further, provided the distribution of radii 
is not too skew, only small errors will be introduced by treating the grain 
radius as having a unique value, a, equal to the mean value of the distri- 
bution. This result suggests furthermore that the shape of the calculated 
A4P—T curve will be insensitive to the shape of the ‘ grains’ considered. 

It is formally possible to make calculations, analogous to those above, 
for metals containing uniformly distributed arrays of line defects (e.g. 
dislocations) or of plane defects (e.g. grain boundaries). In each case, an 
expression similar to (8) is found except that (i) the expressions for K and 
AP’ are slightly modified to take account of the different geometry and (ii) 
in the summand, n is replaced by 3, for line defects and by (2n—1) for 
plane defects, 6, being the ratio of the nth to the first zero of J,(x). 
Further, expressions can be found for the variation of 7’* with « so that, 
provided an experimental value for this variation is known, Q/R and K can 
be calculated and thence 4P/4P’ plotted as a function of 7. An example 
is given in fig. 4 of the shapes predicted by the three different assumptions 
about the nature of the defects for the same set of experimental results. 
Clearly, on changing from point to line to plane defects the width of the 
peak is decreased but the general shape of the curve is not greatly affected. 
In no case has it been possible to distinguish definitely between the three 
theoretical curves on the basis of agreement with experiment. However, 
on physical grounds, while it seems reasonable to consider an oviginally 
uniform array of point defects moving independently, arrays of line or 
plane defects will form networks and the movements of the various parts 
will be interdependent. For this reason, the equations for line and plane 
defects are not set out in detail here. 


§ 3. APPLICATION TO EXPERIMENTAL RESULTS 


In this section, the experimental results of Clarebrough eé al. (to be 
published) are used to determine a series of corresponding values of « and 
T*, From these, Q and D,/a? are evaluated and 4P—T curves plotted. 
The shapes of the calculated curves are then compared with the shapes 
determined experimentally and good agreement is found in the case of the 
low-temperature peaks in nickel. 


Heating rate, « ‘Temperature of first Maximum of APexpt 


(°K/sec) peak, 7'* (°K) (mw) 
0-10 532 61 
0-067 522 42 
0-033 508 25 


The table lists the values of «, 7*, and maximum 4P that are used in the 
calculations, the data being taken from the experimental curves for a set 
of nickel specimens cut from one bar that had been deformed in torsion 
(see fig. 2). When these values are substituted into eqn. (13) and the 
constants adjusted by the methods of least squares and i ageneve approxi- 
mations (as described in the paragraph below eqn. ( 13)), the results are 
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found to fit a straight line, as shown in fig. 3, the final values for the 


- constants being 
Q= 22-6 keal/mole=0-98 ev/defect 


and D,/a*=5-467 x 105 sec}. 
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Solid lines show the low-temperature parts of experimental 4P—7 curves for 
nickel deformed as described in the legend to fig. 1. Heating rates were 
(a) 0-10°c/see, (b) 0-07°c/sec, and (c) 0-03°c/sec. The dashed lines are the 
calculated curves for these cases relative to the dotted horizontal line 
which represents the backward extrapolation of part B of the 
experimental curve. 


Since there are only three points to which the line is fitted, no great 
accuracy can be expected for these figures. However, using the deviations 
from linearity to provide an estimate of the errors involved, Q is found to 
have a standard deviation of +1-0 kcal/mole while D,/a? is accurate to 
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within a factor of 2-8, the large variation in this latter value arising because 
eqn. (13) involves only In (D,/a2). The constants are quoted here to more 
figures than are justifiable since these values were used in the calculations 
to maintain consistency. 

At this stage it should be noted that, for this value of Q, the value of w 
at the maximum of 4P/AP’ is approximately 20 so that the use of the 
asymptotic series (9) and (12) is justified. 

By use of (8) and (10), 4P/4P’ can now be plotted as a function of 7' but, 
in order to compare the results with experiment, we must determine a 
value for 4P’ by comparing the maximum values of AP ip, and 
(AP/AP’)...- There are two possibilities here depending on whether the 
whole or only a part of the experimental rate of evolution of energy at 
the temperatures involved is ascribed to the movement of the point defects. 


Fig. 3 
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Tf the total rate of evolution at 7'* is used, then the calculated AP curves 
are much narrower than the experimental curves. Therefore, it will be 
assumed that the almost ‘ horizontal’ part, B, of the experimental curve 
(see fig. 1) can be produced back to lower bern pereuur and that only the 
energy represented by the area, a’, lying above tus line should be ascribed 
to the movement of the point defects. Use of the corresponding values of 
the maximum of AP,,,, (see table) leads to slightly CE ville for 
AP’ in each case. Therefore, a weighted mean value for 4P’ was she 
mined using weights of 3, 2, and 1 respectively for the runs a sane in ee 
table. This mean value corresponds to Eey=1-22 cal/em*®. The dotte 

curves in fig. 2 show the final plots of era ee the agreement 

d experiment is seen to be reasonapie. 

Bl hi Cian there is some uncertainty about es pee to 
use for 7',, but a consideration of the experimental work stiggents ; ie it 
lies between room temperature and 70°o, this latter temperature being 
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that at which the first evolution of energy occurs. In the curves plotted 
in fig. 2, 7’) has been taken as equal to 70°c (343°K). The choice of 
T)=273°K would have reduced all the 4P values at low temperatures but 


since the decreases range from 25% at 80°c to less than 1% at 120°, © 


the effect on the general shape of the curve is negligible. 
Experimental data for the change in temperature for a peak with 
heating rate is also available for two other cases viz. the high-temperature 
peak in the curves for nickel deformed in torsion (part C in fig. 1) and the 
single peak in the curves for copper deformed in tension. From this data, 
values of Q and K have been calculated and curves for 4P..\. plotted on 
the assumptions of point, line, and plane defects. In fig. 4 a typical set 
of calculated and experimental curves are plotted and it is clear that the 


Fig. 4 
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High-temperature part of the curve in fig. 1 in comparison with curves (dashed) 
calculated on the assumptions of (a) point, (b) line, and (c) plane defects. 
Curve (b) has not been computed for small values of 7 but it must lie 
between curves (a) and (c). 


experimental curve has a much narrower peak than any of the calculated 
curves, a result which is true for both the sets of peaks discussed in this 
paragraph. This lack of agreement is not altogether surprising since both 
these sets of peaks are associated with recrystallization and the present 
model of defects originally distributed uniformly and moving independently 
is clearly not applicable to a metal during recrystallization. 


§ 4. Discussion 


In this section, only the peaks occurring at the lower temperature in 
the curves for nickel will be considered and it will be assumed that they 
are due to the movement of point defects. The nature of these point 
defects has been immaterial up to this stage but the assumption will now 
be made that they are vacancies. Later on, the modifications consequent 


on the alternative assumptions that they are interstitials or vacancy 
pairs will be pointed out. 


— 
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In order to calculate the energy, H, released per defect, an estimate 
of the initial concentration, cy, is required. There is a marked increase 
in density at about the same temperature as 7’* and it seems logical to 
associate the two phenomena. We therefore ascribe the whole density 
change when a sample of the deformed metal is annealed at, say 
(7'*+-50)°k, to the disappearance of the vacancies present in the deformed 
material. The measurements of Clarebrough et al. (to be published) then 
give c)>=2-0 x 10 vacancies/em* from which we deduce that 


H=37 kcal/mole=1-6 ev/vacancy. 


The total energy H+Q=60 kcal/mole should be of the order of the 
activation energy for self-diffusion and although this has not been 
measured directly, Le Claire (1953) gives two empirical methods of 
estimating it, viz. 39x (melting temperature in °K) and 0-66 x (latent 
heat of sublimation). For nickel, these estimates give 67 and 54 keal/mole. 
This excellent agreement may be fortuitous since; if we use Le Claire’s 
theoretical arguments (based on the assumption that diffusion in a 
face-centred cubic metal occurs by the movement of vacancies) to com- 
pute the components of the activation energy of self-diffusion, these 
do not agree nearly so well with the values of H and @. Using Le Claire’s 
formulae, the activation energy for formation of a vacancy in nickel, 
AH,, should be 17-3 kcal/mole and the activation energy for movement, 
AH,, should be 54:8 kcal/mole. These discrepancies will be discussed 
below. 

Le Claire’s work also gives us a means of estimating the order of 
magnitude of Dy. It should be noted that, in the notation of the present 
paper, D, is the frequency factor for diffusion of vacancies already present 
in the lattice and not the frequency factor for self-diffusion by a vacancy 
mechanism since this latter quantity involves the entropy of formation 
of the vacancies. Then, from Le Claire, we have 


D,=vb? exp (4S8/R) and 4S=—Q(u'—p'), . » « « (14) 


where v=frequency of vibration of an atom next to a vacancy, b=lattice 
parameter, 4S=entropy of activation for the migration of a vacancy, 
and p’, p’ are the temperature coefficients of an appropriate elastic 
modulus and the density at the absolute zero of temperature. Using 
(14), with v equal to half the Debye frequency and Young’s modulus as 
the appropriate elastic modulus, Dy can be calculated and hence the value 
of ‘ grain’ radius determined as a=3-4x 10-4 cm. a Eee 

In view of the possible errors in the value of D)/a? and in view of the 
assumptions made in deriving and applying eqn. (14), it seems wbwase 
to attach too much weight to this value of a. However, this oe 
size’ of 7 is more comparable with the limits of ‘ particle dimes Viz. 
1-9 to 0-041, found by Gay and Kelly (1953) in nickel rolled 33° than 


with the crystal size of approximately 100p that Clarebrough et al. have 
ion of their specimens. On this basis 


» found by microscopic examinatl , 2 3 ba 
the * fae ’ of the present theory will be tentatively identified with the 
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‘particles’ of Gay and Kelly and it will therefore be assumed that 
most of the vacancies disappear by absorption onto the edge dislocations 
which exist in the ‘ particle’ boundaries (Gay, Hirsch and Kelly 1954). 

On this interpretation of the results, H becomes a measure of the mean 
energy evolved when a vacancy is absorbed onto an edge dislocation in 
a particle boundary and this should not be identical with the 4H, of 
Le Claire’s theory since this latter is the energy for formation of a vacancy 
in a perfect lattice. A discussion of the relative magnitudes of H and 
AH, is pointless since the value of H depends directly on the value 
assumed for cy i.e. on the assumptions made about the relationship 
between overall density change and number of vacancies originally 
present. 

The discrepancy between @ and the theoretical estimate of 4H, is 
more serious since both purport to be measures of the same quantity, viz. 
the activation energy for migration of a vacancy in nickel. This is not 
strictly true since @ has been determined from measurements on a 
deformed metal and the work of Gay et al. has shown that even within 
the ‘ particles ’ there is still a certain amount of strain. Thus, we might 
expect @ to be less than the value 4H, found from work on annealed 
metals. Furthermore, Le Claire’s formula for 4H, contains two possible 
sources of error viz. (i) a numerical factor calculated from diffusion data 
on other metals and (ii) an uncertainty in the choice of an appropriate 
elastic modulus. In this connection, it is interesting to note that Le 
- Claire’s formula gives 4H, for copper equal to 36 kcal/mole. This is 
considerably higher than the theoretical value of about 22 kcal/mole 
given by Huntington and Seitz (1942) or Bartlett and Dienes (1953). 

Let us now consider the possibilities that the moving defects are either 
interstitials or vacancy pairs. There are no theoretical values for the 
activation energies of migration for either of these defects in nickel but 
for copper, Huntington (1953) estimates 0-25 ev for the energy of migra- 
tion of an interstitial atom while Bartlett and Dienes (1953) estimate 
0-34 ev for the energy of migration of a vacancy pair, both values being 
considerably less than @. While Huntington’s work cannot easily be 
extended to nickel, substitution of the appropriate figures for nickel into 
the formulae of Bartlett and Dienes causes very little change from the 
copper results. Thus the measured value of Q provides strong support 
for the suggestion that the defects are vacancies. 

A further argument against the postulate of interstitials causing the 
peak is provided by the measurements of density. On the same assump- 
tion as before, i.e. that the first increase in density is associated with the 
peak in the 4P curve, this implies that the disappearance of an inter- 
stitial atom increases the density. While this seems unlikely at first 
sight, it is not impossible since the removal of an interstitial atom does 
allow the lattice to become more close-packed.’ However, since atoms 
are not hard spheres, the elastic forces will almost certainly ensure that 
any density increase due to the removal of an interstitial atom is less 


Release of Energy from Deformed Nickel 97 


than 10% of that associated with the removal of one vacancy. This 
would imply #<0-16 ev per interstitial, a value that seems far too low 
to be the energy released on the disappearance of an interstitial atom. 
In the case of vacancy pairs, the density argument follows the same 
lines as for single vacancies and we find #=3-2 ev per vacancy pair which 
is not unreasonable. 

Therefore, it can be concluded that the assumptions set out in § 2 
provide a model that gives a consistent description of part of the release 
of energy from deformed nickel. Furthermore, while it is not impossible 
that the defects involved are interstitials or vacancy pairs, the data 
deduced from the experiments can be interpreted most reasonably on 
the assumption that the defects are single vacancies. 
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XIV. CORRESPONDENCE 


The Vibrations of a Stressed Framework 


By H. C. Lonever-Hicens 
The University Chemical Laboratory, Cambridge 


[Received October 1, 1954] 


In this note equations are derived for the normal modes of vibration of a 
system of massive particles connected in pairs by light springs. The case 
of interest is that in which the number of springs exceeds the number of 
internal degrees of freedom of the system ; for then in equilibrium the 
springs will be under stress, and the vibrational modes will depend on the 
equilibrium stresses as well as on the force constants of the springs. 

Let r, be the position vector of particle « in equilibrium, and denote 
r,—t, by r,3- Then the force exerted on particle « by particle f in equili- 
brium is 

— dup Lem 
ep 
where ¢,¢ is the equilibrium tension in the spring «f. 

Now let the system be distorted slightly from equilibrium, so that the 
position of « becomes r,+s,. The vector distance of « from f is now 
Tp +Sup, Where $S,,—S,—S,;, and the scalar distance is r,,;+ 9,3, where 
Tah Iap=Tep + Sug The tension in the spring «6 may therefore be written 
dup t Yup Iaps Where yj, 1s the force constant, and the force exerted on « by B 
is 

=z (Pap t Pup Yap) (Tap t+ Sag) 
| Top Wap 
which, to the first order of small quantities, is equal to 
Ww 
{ata meget T Xap sa} 


where Xap= bap] "aps Wap=Vup— Pap] Tap: 
Now in a normal mode of frequency v, the total force on « must equal 


: Fa 
—4n*v*m,S,, where m, denotes the mass of the particle. But 2 y,, h3—O 


since the total force vanishes in equilibrium. Hence, if A=47*22, 


cogs ( r, 
AM, Se et {ewe ha + Xap so} : 
J B Yop 
Therefore 


AS,,= X(S,—S$,,) 


1 FO Wop Vup Vp 1 ao Ww, r. 
——— Dk {caters + Xap Su) m, Pe ee dye a eve s¢ 
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Taking the scalar product of both sides with Txy/Tay, One obtains 
1 #2 
vy =a { Poy » Szp Lee Levee 3 
ay Me ; Wap yap dap Xap pa + Mm, aa Wye Cove Guat Xue : oe x ’ 
where ¢,,, is the cosine of the angle ya8. Now 


Pay - Syg=(Tap +8 py) - Sag 
=T1ep Wap thay « (Sxy+S,g) 
= Io6 Qua + (Vpa t Tap) Say—"By Upy 
=p Tap —" By Ipy Tay Tay Tap + Sy: 


” } Fa 
Using again the fact that 2’ y,, r,,=0, we obtain finally 
B 


] #4 


— Tes r 
A My 4 fo Crop Yap + Xap (“2 Vap— “8 +4) 
1 * 


a é P coe 
a m, 2 {. Cope Dy Xye ( Tye 7, at dy )} : 
This is a set of homogeneous linear equations of the form 
a Ais G=W 
a 


where i and j denote pairs of particles. Ifi andj refer to the same pair of 
particles («, y), then 


1 Fa il #Y 
A= (yt E nas) + a (Hort 2 x) 
If i denotes the pair («, y) and j the pair («, 8) then 
1 Top il "ap 
A= Mr (ox, Coayt XaB =) ak m, XBy Pay’ 


and if i and j have no common member then A,;=0. 

If there are n particles, there will be $n(n—1) equations altogether, 
determining 3n—6 modes of vibration, if the system is non-planar. The 
difference between these two numbers is the number of geometrical 
relations between the q;. For there to be permanent stresses in the 
system this number must be positive, which implies n>4. However, if 
the equilibrium configuration is planar (or linear) the equations give only 
the in-plane (or in-line) vibrations. This is because the out-of-plane or 
out-of-line vibrations do not-alter the distances between particles to the 
first order. In the absence of permanent stresses these vibrations would 
consequently have zero frequency ; but in the presence of stresses this is 
not so, as the following argument shows. 

Suppose the equilibrium configuration is planar, and cons 
in which every particle moves a small distance perpendicul 
The extensions q,,3 vanish to the first order, and so the force on « reduces to 

Fa 
i Xees Sass 


ider a distortion 
ar to the plane. 
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where S,, is perpendicular to the equilibrium plane. Therefore 
Fu ee 
Say= Mx 2 Xap Sap— m, = Xye Sye 


taking the scalar product of both sides with s,,/s,,, one obtains 
] #a« ] = 
Say= Ma a Xap Sop Mm, a Xy Sye 


which is a set of equations of the form 
2 Bj; 8;=A8;;, 
j 


like the equations for the in-plane modes. A necessary condition for the 
out-of-plane frequencies to be finite is therefore that the permanent stresses 
shall not all. vanish. ; 


. Angular Correlation in the Reaction *N (d, ayy) ?C 


By J. SEED 


Cavendish Laboratory, Cambridge 


[Received November 8, 1954] 


AN alpha-particle group in the reaction #4N (d, «) 12C, corresponding to a 
level in carbon-12 at 7-5 Mev, was first reported by Holloway and Moore 
(1940) ; recently the existence of the alpha-group has been confirmed by 
the more accurate work of Dunbar et al. (1953), who assign an energy of 
7-68 mev to the 2C level. 

Evidence for the existence of this level has also accumulated from 
studies of the reaction *Be («,n) 12C. Using polonium alpha particles, 
Harries and Davies (1952) observed electron—positron pairs of total energy 
7 Mev in a cloud chamber. No gamma rays of energy around 7 Mev have 
been observed, although Beghian et al. (1953) and Uebergang (1954) have 
reported 3 Mev gamma rays in coincidence with the 4:5 Mev gamma rays 
from the first excited state of carbon-12. The above evidence alone 
might lead one to conclude tentatively that the 7-68 mev level has spin 
either zero or >4. 

In the present author’s experiments, the angular correlation between 
3Mev and 4:5Mev gamma rays has been measured in the reaction 
AN (d, ayy) *C. Targets of tantalum nitride, made by heating tantalum 
metal in dry ammonia, were bombarded by a 1 pA beam of 630 kev 
deuterons. Gamma rays emitted by the target were detected by two 
sodium iodide scintillation counters, the angle between which could be 
varied in the plane perpendicular to the incident deuterons. The output 
of the one counter was fed into a single-channel kicksorter adjusted to 
select pulses corresponding to the energy of 4-5 Mev gamma rays less two 
annihilation quanta ; the output of the other counter was fed through a 
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gate into a 60-channel pulse height analyser (Hutchinson and Scarrott 
(1951)). The gate operated when a fast coincidence from the two 
y-counters was accompanied by a pulse from the output of the single- 
channel kicksorter. The coincidence resolving time was 3-5 10-8 sec 
using a circuit of the type described by Bell, Graham and Petch (1952). 

In the experiments, coincidence y-ray spectra were observed on the 
60-channel pulse analyser for five different angles between the gamma- 
counters. At each angle, peaks at approximately 2-5 mev and 2 Mev 
appeared, corresponding to the one-escape and two-escape annihilation 
peaks of a 3 Mev gamma ray in coincidence with the 4:5 Mev gamma ray. 
These peaks were superimposed on a background chiefly due to 1-5 Mev 
and 2 Mev gamma ray coincidences with the 5 Mev gamma rays in the 
reaction 14N (d, pyy) °N. For the purposes of the correlation, the in- 
tensity of the 3 Mev gamma ray was estimated at each angle from the 
strength of the ‘ one-escape ’’ peak, monitored on the number of gating 
pulses applied to the 60-channel analyser. At all angles this number was 
much greater than the number of ‘ one-escape ’ pulses, owing to the pre- 
sence of the background coincidences from N (d, pyy) /N. 

In order to check that the 2-5 Mev peak arose solely from coincidences 
with the 4-5 Mev gamma ray, a second experiment was performed in which 
the single channel kicksorter in series with the one counter was adjusted 
to detect only gamma rays of energy greater than 4-5 Mev. No pulses at 
2-5 Mev were observed in this experiment, although a small peak at 2 Mev 
was still evident, corresponding to the total capture energy of'a 2 Mev 
gamma ray in 1N (d, pyy) ®N. 

The results of the coincidence measurements (about 200 hours’ running 
time) are shown in fig. 1. The full curve represents 

I(0)=1—8 cos? 0+-4 cos* 8 

corrected for the finite solid angles of the counters, and it is seen that 
within the limits of accuracy of the experiment the curve is a reasonable 
fit to the experimental points. This curve would result if the 7 ‘68 Mev 
level had spin and parity (0+), giving a pure 0-- > 2-- 7 oF transition 
involving two electric quadrupole gamma rays. Odd parity for the level 
would give a similar result, however, and higher spins such as 3-++ or 4+. 
cannot be definitely ruled out. f : 

It has been pointed out (Salpeter 1953, Opie 1952) that, in stars that 
have largely exhausted their central hydrogen, oxygen and carbon can be 
produced via the reactions 

4He+‘4He — ®Be 
8Be-+4He — ?C 
122C 1 4He — 16O 


Hoyle (1954) has concluded from the observed cosmic abundance ratios 
of the elements that the latter reaction should havea pon ene at 0-31 Mev, 
corresponding to a level in carbon-12 at 7-68 MeN: If he devel has edd 
parity with zero spin, however, the reaction §Be--4He — 3°C cannot take 


place. 
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Fig. 1 
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An Estimate of the Mean Lifetime of the 0°-particles 


By D. I. Paar* 
The Physical Laboratories, The University, Manchester} 


[Received October 15, 1954] 


AN estimate of the mean lifetime of the 6°-particles has been made by a 
method similar to that previously used for the A°-particles (Page 1954). 
From a sample of 357 V°-decays obtained at the Jungfraujoch (3 580 m), 
14 events were selected which satisfied the following criteria : 

(a) The momentum of the positive secondary particle and the ionization 
density were consistent with the particle being a light meson and incon- 
sistent with its being a proton, 

(5) the dynamics of the decay were consistent with the scheme 
(Thompson 1954) 

0° > at++a—+214 mev. 


It is estimated that less than 10°/ of these events are due to the decay of 
particles other than @°-particles (Discussion at the Bagnéres de Bigorre 
Conference, 1953). In all cases the potential path length of each decaying 
particle was limited by the length of the positive secondary track which 
must be allowed for the momentum of the particle to be determined. 
This in turn was limited by the maximum detectable momentum measure- 
able in the chamber. 

An analysis of the actual and potential times of flight of the particles 
(Bartlett 1953) gives the lifetime as 


+= (0-740:3) x 10-1 sec. 


Combining this value with all the previously published results (Gayther 
1954) the weighted mean value for the lifetime becomes 


7=(1-579:4) x 10-1 sec. 


I am grateful to the administration of the Hochalpine Forschungs- 
station, Jungfraujoch, for providing research facilities at the Jungfraujoch 
and to the Department of Scientific and Industrial Research for a main- 
tenance grant. 
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Reversible Generation of Temperature Waves (‘ Second Sound ’) 
in Liquid Helium II 


By N. Kurri and J. McIntosH 
Clarendon Laboratory, Oxford 


' [Received December 10, 1954] 


THE existence of temperature waves in liquid helium IT was first postu- 
lated by Tisza (1938). Landau (1941), on the basis of his quantum 
hydrodynamical treatment of liquid helium, proposed the existence of a 
phenomenon which he called the ‘second sound’. An unsuccessful 
attempt to detect this effect by an acoustic method, namely by the beats 
in standing waves radiated by an oscillating piezo electric quartz crystal 
was made by Shalnikov and Sokolov (quoted by Peshkov 1947). An 
analysis of this negative result led Lifshitz (1944) to predict the essentially 
thermal nature of ‘second sound’. Peshkoy (1944), using thermal 
generation succeeded in demonstrating the existence of ‘ second sound ’ 
which is identical with the temperature waves proposed by Tisza. Since 
then many experiments have been done on these temperature waves, 
but with the exception of Peshkov’s (1948) method of ‘ filtration of the 
superfluid ’ they appear to have been always generated by irreversible 
Joule heating. In this method there is a net heat flow and heat dissipa- 
tion in the liquid helium and this might affect the results especially at 
the lowest temperatures. 

The temperature oscillations of a suitable paramagnetic salt placed 
in an alternating magnetic field provide a method ‘ par excellence ’ for 
producing second sound reversibly.* This method should compare 
favourably with the filtration of the superfluid method which becomes 
ineffective at low temperatures. We have therefore carried out some 
experiments between 1-3°K and the A-point in order to prove the feasi- 
bility of generating second sound by means of a paramagnetic salt. 

Standing waves of second sound were produced in a resonant cavity 
of 2 em diameter and 2-3 em length. A carbon resistance thermometer 
which served as detector formed one end of the cavity while the other 
end was formed by the smooth end-face of a cylindrical pill of compressed 
ferric ammonium alum powder. An alternating magnetizing field of 
up to 120 gauss at frequencies between 150 and 1200 cycles per second was 
produced by a small external coil cooled in liquid hydrogen. The un- 
wanted signal produced in the thermometer circuit by the stray field of 
this coil was compensated by injecting in the thermometer circuit a 
voltage of suitable phase and amplitude. To obtain maximum sensitivity 
and to eliminate the effect of any residual pickup, the output of the 


* This method was probably first considered in the autumn of 1938 when 
possible experiments to demonstrate his temperature waves were discussed 
with Professor Tisza. 
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amplifier, normally displayed on a cathode ray oscilloscope, could be 
coupled to a phase sensitive detector and the amplitude of the signal 
measured on a meter. 

The results may be summarized as follows. Standing waves were 
_ observed with the fundamental frequency and up to the third harmonic, 
By measuring the resonance frequency at selected temperatures the 
temperature dependence of the velocity of ‘second sound’ was also 
determined. The values found agreed within the limit of experimental 
accuracy with those of previous authors. The temperature excursions as 
measured on the carbon resistance thermometer were only of the order 
of 10-*°, although the temperature amplitudes of the spin system of the 
paramagnetic salt were of the order of 10-2°. This is due partly to the 
effect of the spin lattice relaxation time but chiefly to the relatively large 
specific heat of the liquid helium in this temperature range. It is worth 
pointing out that with a Joule heater a power input of the order of 10-2 
watts would be required to obtain temperature amplitudes of 10-4° in . 
liquid helium. 

Experiments were also carried out in a superimposed constant magnetic 
field of 600 gauss so as to increase the temperature amplitude for a given 
a.c. field. This method has, however, the slight disadvantage that the 
frequency of the second sound is now the same as that of the a.c. field 
which makes it more difficult to distinguish spurious effects from genuine 
resonance. With no superimposed steady field the second sound fre- 
quency is twice that of the a.c. field frequency. 

Our experiments have shown that the efficacy of this method of 
generating ‘second sound’ increases with decreasing temperature. In 
other words, any decrease in amplitude due to increasing spin lattice 
relaxation time is outweighed by the decrease of specific heat of liquid 
helium. Moreover it should be possible to find substances with shorter 
spin lattice relaxation times than that of ferric ammonium alum. It 
seems therefore that this method will be suitable for the temperature 
range below 1°k. Experiments with only very small heat dissipation 
might well help towards a better understanding of the exact nature of 
‘second sound’ in this temperature range. The paramagnetic method 
could also be used for the pulse technique and it will thus be possible 
to generate ‘ cold’ pulses as well as ‘ heat.’ pulses. 


This work was done during J. McIntosh’s tenure ot a Berry Scholarship 
of the University of St. Andrews for which grateful acknowledgement 
is made. 

REFERENCES 
Lanpav, L., 1941, J. Phys., U.S.S.R., 5, 71. 
Lirsurrz, E., 1944, J. Phys., U.S.S.R., 8, 110. i 
Prsukov, V. P., 1944, J. Phys., U.S.S.R., 8, 381; 1947, Rep. Intern. Conf. 
Cambridge, Phys. Soc., Vol. II, 19 ; 1948, J. Lup. Theor. Phys., USS it., 


18, 867. 
Tisza, L., 1938, Comptes Rendus, Paris, 207, 1035, 1186. 


106 Correspondence 


Spin Orbit Coupling and the Mesonic Lamb Shift 


By E. W. Larne 
Department of Natural Philosophy, University of Glasgow 


[Received November 18, 1954] 


Ir has recently been shown by Chisholm and Touschek (1953) that the 
self-energy corrections for a nucleon moving in a scalar potential well 
U(r) lead to a strong spin orbit coupling for pseudoscalar mesons with 
pseudoscalar coupling. Assuming charge symmetric meson theory, 
these authors obtain a spin orbit coupling 


AU*®-=3(g/40Mf)? = Gg. L, Corp. le Cerca) 
This is of the right order of magnitude, but of the wrong sign, for applica- 
tion to the nuclear shell model. 

It is the purpose of this note to show that a more detailed investigation 
of the total matrix element reveals contributions to the spin orbit coupling 
from the small components of the initial and final spinor wave functions. 
These contributions are such that they reverse the sign of the spin orbit 
coupling, while leaving its magnitude unchanged. 

We shall employ the usual Feynman-Dyson prescription for evaluating 
the S-matrix element A(p., p,;) in momentum space. The only second 
order graph which contributes is shown in the figure. 


k 
Tie Sa ae 
ge re 
i \ 
/ . 
———1__><_>—__1__5 — 
P, p,tk U@) pot k p> 


This leads to the following expression for the matrix element 


Alo, | - = : = 
(Po Pr (P+ p2)[(p, Lk MI (p,-- EP] 


| (2) 
where Ultr)= | ge Utes, 
A(P», Py) is logarithmically divergent. By the usual renormalization 


procedure, A)= A(po, po) is subtracted from A(p5, p1). (Po in Ay is a 
free particle momentum.) 


js Sige ap, PDadvalty-(Patk)— MIU Qlty.(prtk)—Mysh(p1) 
Crayne 
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The finite part A,= A— A, contains four terms which contribute to 
the spin orbit coupling, namely : 


A,=A, (2) Y iV sP2iP1z ¥(P1) U(Q), 

A,=A, ii po) vaVi Po; (D1) UY), (3) 
A,=A, (pe) y;P3; u(p1) U(Q); 

A,=A, (Do) u(p) U(q), 


where 7,j=1, 2, 3 and summation is implied. Here P ni IS a linear 
combination of p,; and p5;, say a,)1;+6;,P0;, k=2, 3. The constants 
A,, Az, Az, Ag, 4,, 6, are easily obtained by standard methods. 

The A,(u=1, 2, 3, 4) contain the following terms, A,’, on reduction 
to two-component spinors y in the non-relativistic limit, 


Ay'=1A,xto . (P,x px) U(q)x | 
=1A,S, say 
A,’=(b,—4) A,S/2M, Mah of eae) 
A;'=(a3+-63) AgS/2M, 
A,’ =—1/(2M)? AS, J 
rahe 4 : 
us ae ak L, in configuration space. 
The results obtained are as follows :— 
1dU 
| facta t Dnata s site ie 
A,'=3(9/47M) ae L, a 


A,’ is negligible, since, to the order considered, calculations show that 
the static radiative correction amounts to only GY, Of U(r). 
Thus, the total contribution to the spin orbit coupling is 


1dU 
—3(g/4n M1)? — o.L. cts pyle pee) 


dr 
This is just a change of sign from the previous work by, Chisholm and 
Touschek, referred to above, which reported the contribution of the 

‘ alone. 

eee clis obtained are encouraging insofar as they afford a possible 
field theoretical picture for the nuclear shell model. One must bear in 
mind, however, that a systematic search for spin orbit coupling from all 
possible sources should be made and that only the total contribution 
should be compared with the experimentally observed spin orbit coupling. 
For example, Klein (1953) obtains a spin orbit coupling of the wrong sign 
for the shell model as a velocity dependent correction to fourth order 
two-body forces. Further work is being done by the author along the 


lines indicated above. 
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In conclusion the author wishes to extend grateful thanks to Professor 
J. C. Gunn and Dr. J. 8S. R. Chisholm, under whose guidance this work 
was performed. 
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ERRATA 


On the Evaluation of Characteristic Temperatures for Cubic Crystals, by 
A. B. Buatta and G. E. Taussr, 1954, Phil. Mag., 45, 1211. 


(1) In eqn. (3) read 47a 9K, instead of 47a. 


(2) Multiply the right-hand side of (5) by (p)3/2, where p is the density 
of the substance. (The @-values given in the table were, however, 
calculated with the correct formula.) 


ie They 


XV. REVIEWS OF BOOKS 


Atlas of Typical Expansion Chamber Photographs. By W. Genrner, H. 


Mater-Lerpyirz and W. Boraz. [Pp. 199.] (Pergamon Press Limited.) 
Price £5 5s. 0d. 


Tuts volume forms a completely new and up-to-date version of the ‘ Atlas 
typischer Nebelkammerbilder ’ published by the same authors in 1940. An 
impressive selection of photographs, many as yet unpublished, illustrates the 
numerous contributions made by cloud chamber observations to the field of 
nuclear physics. The subjects covered include natural and artificial radio- 
activity, nuclear transmutations, and the properties of the various types of 
fundamental particles and radiations. The atlas is completed by a small but 
representative collection of cosmic-ray events—these having already been 
comprehensively surveyed in “Cloud Chamber Photographs of the Cosmic 
Radiation ’ by Rochester and Wilson. 

The quality of reproduction of the photographs is excellent. They are 
adequately explained by legends in English, French and German, and supple- 
mented by some thirty useful diagrams and tables. Not only will this volume 
be invaluable to the student of nuclear physics, but it will appeal also to the 
non-specialist, for the photographs display in a straight-forward and vivid 
manner the fundamental processes resulting from interactions in the atomic 
nucleus. ETE 


Rocket Exploration of the Upper Atmosphere. Edited by Dr. R. L. F. Boyp 
and Dr. M. J. Seaton in consultation with Prof. H. 8. W. Massry, F.R.S. 
(London: Pergamon Press Ltd.) [Pp. viii+376.] Price 75s. 


Workers in America have opened up a new method of attack on the problems 
of upper atmospheric physics by the use of rockets, of the type devised by the 
Germans during the war of 1939-45. The papers presented in this volume 
constitute a good survey of the results obtained up to August 1953. The 
greatest height reached seems to have been about 400 km, while a height of 
150 km was reached several times. 

There is always some doubt about deductions made in a very tenuous 
atmosphere from a rocket which inevitably carries some gas with it, and which 
is moving at a velocity comparable with the gas kinetic velocity of the molecules, 
and some of the deductions are still open to doubt. In spite of that doubt 
the results have already made a considerable difference to the theory of the 
ionosphere. The most direct measurements are those of the intensity of 
solar x-radiation at different heights. Less direct, but very important, are 
those in which the density and temperature of the atmosphere are measured. 
The results of these experiments have already revolutionized our thinking about 
the radiations responsible for ionizing the E layer, and about the air density 
in the F region. 

In studies of the composition of the upper atmosphere the techniques have 
‘included direct sampling, and the use of a mass spectrograph in the rocket. 
Direct measurements of electron and ion density in the ionosphere have been 
attempted. All these experiments are very difficult and we cannot yet be 
certain that the results are unambiguous. ; 

Tn one of the most direct, and important, experiments the earth’s magnetic 
field was measured at different heights and the results seem to show that the 
currents of the ‘ atmospheric dynamo’ theory of geomaguetism flow at a 

i 0 km. 
ac ean be of great value to all students of the upper nals oe 
particularly because much of the American work has up to now only appearec 


4 } Rape BE: 
in progress reports and has not found its way into acknowledged hi ~ ‘ hes 
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Quantum Mechanics. By F. Manpu. (Butterworths Scientific Publications.) — 


[Pp. 232.] Price 35s. 


Non-RELATIVISTIC quantum mechanics seems here to stay, and most physicists 
must understand it in the same way that engineers understand dynamics ; 
that is, they must know how to calculate with it to solve specific problems. 
Dr. Mandl’s treatment is in this spirit. All the usual topics are included 


(except field theory), but only up to the degree of complexity necessary to — 


demonstrate an important principle or technique: the reader can work the 
heavier algebra for himself, in numerous examples with detailed hints. The 
author writes: “‘ Apart from classical physics, calculus, and vector analysis, 
the reader is assumed to be familiar with the qualitative concepts of quantum 
theory, such as the wave nature of matter, and their experimental basis ”’. 
The special virtue of the book is that without attempting real mathematical 
and metaphysical rigour (although there is an illuminating note on the Dirac 
6 function), it sets forth clearly and concisely the various viewpoints, formula- 
tions, and notations of such men as Schrodinger, Heisenberg, von Neumann 
and Dirac in a logical whole, so that the student will have knowledge and 
experience of all the different techniques when he comes to read and to write 
original research papers. This eclectic, up-to-date account (it even shows how 
group theory stems naturally from atomic physics) should make a most useful 
textbook, suitable for good students in their third year and for all post-graduate 
students, whether they specialize on experimental or theoretical physics. 


Transactions of Symposia on Applied Mathematics, Vol. I. [Pp. 243.] (Inter- 
science Publishers.) Price $5. 


Tx1s volume contains a collection of papers delivered at a symposium on 
fluid mechanics and computing held at New York University in April 1953 and 
sponsored by the American Mathematical Society and the U.S. Army Office of 
Ordnance Research. 

In all there are 14 papers by various well-known scientists. They consist 
of a paper on the theory of homogeneous turbulence, by G. Birkhoff, a critical 
discussion of approximation techniques used in boundary layer problems, by 
G. F. Carrier, four papers on computational methods by J. H. Giese, M. Lotkin, 
P. D. Lax and L. H. Thomas, and eight papers on compressible flow by M. J. 
Lighthill, G. 8.8. Ludford and M. H. Martin, L. Bers, A. Weinstein, P. Germain, 
R. von Mises, G. E. Hudson and H. Schardin. 

Of particular interest are the four papers on the problem of transonic flow, 
which give a good picture of the present state of knowledge on this difficult 
problem, and contain some interesting speculations. We, 


BOOK NOTICES 
Probability and Information Theory with Applications to Radar. By P. M. 
Woopwarp. (London: Pergamon Press Ltd.) [Pp.128.] Price 21s. 


The Atmospheric Lunar Tides. By Ryuxicrt Sawapa. (New York: Univer- 
sity Press.) Price $1.50. 

British Journal for the Philosophy of Science. Vol. V, No. 19. (Edinburgh : 
Thomas Nelson and Sons Ltd.) Price 7s. 6d. 


Nuclear Species. By H. E. Hunrtey. (London: Macmillan & Co. Ltd.) 
[Pp..193.] Price 21s. 


[T'he Hditors do not hold themselves responsible for the views 
expressed by their correspondents. | 


